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METHOD AND SYSTEM FOR PREDICTING NUCLEIC 
ACID HYBRIDIZATION THERMODYNAMICS AND 
COMPUTER-READABLE STORAGE MEDIUM FOR USE THEREIN 

BACKGROUND OF THE INVENTION 

5 1 . Field of the Invention 

This invention relates to methods and systems for predicting nucleic 
acid hybridization thermodynamics and computer-readable storage medium for use 
therein. 

2. Background Art 

10 Improvement of the efficiency of hybridization-based techniques 

requires the optimization of the binding between two sequences. Accurate 
prediction of the thermodynamics allows optimal choice of the sequences, 
temperature, and salt conditions. Hence, the prediction of nucleic acid 
thermodynamics is important to optimize techniques like PCR (Saiki et aL, 1988), 

15 Southern and Northern blotting (Southern, 1975), antigene targeting (Freier, 1993), 
and Kunkel site-directed mutagenesis (Kunkel et aL, 1987). 
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Hybridization prediction is also important for designing DNA 
microchips that have a wide field of application ranging from diagnostics (Hacia, 
1999; Yershov et al., 1996) to gene expression analysis (Ferea et al., 1999) and 
drug discovery (Debouk and Goodfellow, 1999). Microchips contain a large 
5 number of DNA probe sequences that have to be designed to specifically hybridize 
target sequences in a pool of DNA fragments. First, a DNA probe should be 
designed to bind to only one site of only one DNA target. Second, the different 
DNA probe sequences need to hybridize to their targets under the same temperature 
and solution conditions. Moreover, in sequencing by hybridization (Fodor et al., 
10 1993; Mirzabekov, 1994) where microchips are used to determine the sequence of 
given DNA, one has to be able to know hybridization thermodynamics to 
discriminate signals resultmg from perfectly matched and mismatched probe/target 
hybridizations. 

Another widely used technique that requires hybridization prediction 
15 is the fluorescence in situ hybridization (FISH) technique (Gall and Pardue, 1969). 
In this technique, a fluorescent tagged nucleic acid probe is designed to specifically 
hybridize cellular or tissue section nucleic acids. The target of these probes can 
either be endogenous DNA, messenger RNA or viral and bacterial sequences. 

Therefore, FISH is used to monitor gene expression (McNicol and 
20 Farquharson, 1997), detect infectious agents (Bashir et al., 1994; McNicol and 
Farquharson, 1997; PoUanen et al., 1993), study cell cycle (McNicol and 
Farquharson, 1997), map chromosomes and study nuclear architecture (Heng et al., 
1997). It was also determmed that a set of probes can be used simultaneously 
(multiFISH) to detect different loci (Pagon, 1997). Once again, prediction of 
25 hybridization is essential to insure specificity. Nucleic-acid hybridization prediction 
is also important for the design of oligonucleotide aptamers or antisense 
oligonucleotides (Cohen, 1992) that can be used for various therapeutic applications. 
A new type of probes known as molecular beacons (Bonnet et al., 1999; Tyagi et 
al. 1998) that are very specific has been developed and shown to be efficient for 
30 mutation analysis (Giensendorf et al., 1998) and multiplex detection of single 
nucleotide variations (Marras et al., 1999). The design and prediction of the 
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thermodynamics of these beacons is helped by hybridization thermodynamics 
prediction (Bonnet et al., 1999). Accurate prediction of hybridization is also 
important for the practical realization of DNA-based or more generally nucleic acid- 
based computers. (Adleman, L.M., 1994), 

5 The development of molecular biology techniques based on 

hybridization (PGR, FISH, DNA microchips, etc.) has resulted in a need for 
efficient automated ways to design probes and primers. In the last decade, 
numerous algorithms have been developed to optimize the design of primers and 
probes for various applications (Rychlik and Rhoads, 1989)(Breslauer et al., 1986; 

10 Chen and Zhu, 1997; Dopazo et al., 1993; Haas et al., 1998; Hillier and Green, 
1991; Hyndman et al., 1996; Li et al., 1997; Link et al., 1997; Pesole et al., 1998; 
Proutski and Holmes, 1996). Numerous unpublished software to predict primers 
are also made available by research groups and biotech companies on the World 
Wide Web (PrimerS from the Whitehead Institute for Biomedical Research. Primer 

15 Express™ from PE Biosystems, DNAstar from IDT, etc.) 

There are currently many software packages on the market for DNA 
primer design including: OLIGO, PRIMER PREMIER, OSP, GCG, PrimerMaster, 
and Primo. None of the current programs, however, were written by experts in 
DNA thermodynamics; thus, there are many improvements that can be made. 

20 Nearly all of the current software packages contain mistakes that result from a lack 
of understanding of the underlying theory of DNA hybridization. PGR is a fairly 
robust process and thus even crude programs make predictions that work 90-95 % 
of the time. Multiplex PGR primer design, however, is not at all trivial and 
detailed knowledge of the physical chemistry of DNA hybridization, and the 

25 availability of an accurate thermodynamic database are essential to reliable design 
of multiplex PGR primers. In multiplex PGR, several primers must be designed to 
specifically bind to different sites on target DNA at a given temperature with 
minimal background binding to mismatch sites and with minimal cross- 
hybridizations between pairs of primers. The design of molecular beacons for DNA 

30 oligonucleotide arrays is also very challenging because of the complex competing 
equilibria. 
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Most of the existing programs aiming at finding an optimum probe 
that binds a specific location on a target, however, do not include accurate stability 
rules for hybridization and neglect or poorly approximate competitive binding sites, 
strand folding and strand dimerization. 

5 U.S. Patent Nos. 5,593,834 and 6,027,884 to Lane et al. disclose 

methods to design and construct DNA sequences with selected reaction attributes. 

In summary, prediction of nucleic acids thermodynamics is important 
to optimize various molecular biology techniques including multiplex PGR, DNA 
microchips, molecular beacons, and fluorescence in situ hybridization. Most of the 
10 available programs for probe design do not include a complete parameterization and 
often do not account for mismatches. Moreover, single strand folding is not taken 
into account, which often leads to inaccurate predictions. 

SUMMARY OF THE INVENTION 

An object of the invention is to provide a method and system for 
15 predicting nucleic acid hybridization thermodynamics and computer-readable storage 
medium for use therein wherein the invention utilizes a thermodynamically rigorous 
approach to evaluate the quality of probes and simulate probe/target hybridization. 

Another object of the invention is to provide a method and system for 
predicting nucleic acid hybridization thermodynamics and computer-readable storage 
20 medium for use therein wherein the invention also takes into account single strand 
folding thermodynamics to calculate effective hybridization thermodynamics. 

In carrying out the above objects and other objects of the present 
invention, a method for predicting nucleic acid hybridization thermodynamics is 
provided. The method includes providhig a database of fhemiodynamic parameters, 
25 receiving hybridization information which represents at least one sequence, 
receiving correction data, receiving a first set of data which represents hybridization 
conditions, and calculating hybridization thermodynamics including net 
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hybridization thermodynamics based on the hybridization information, the 
thermodynamic parameters, the correction data and the first set of data. 

The hybridization thermodynamics of individual single stranded, 
bimolecular and higher order complexes may be statistically weighted in a numerical 
5 process and the equilibrium concentration of each species is output. 

The correction data may include folding correction data and/or linear 
correction data. 

The thermodynamic parameters may include DNA thermodynamic 

parameters. 

10 The DNA thermodynamic parameters may include dangling end 

parameters and/or coaxial stacking parameters. 

The DNA thermodynamic parameters may further include terminal 
mismatch parameters . 

The thermodynamic parameters may include RNA thermodynamic 
15 parameters and/or hybrid DNA/RNA thermodynamic parameters. 

The thermodynamic parameters may further include DNA loop 
thermodynamic parameters. 

The hybridization information may represent top and bottom strand 
sequences which form a duplex and wherein the hybridization thermodynamics are 
20 calculated for the duplex. 

The hybridization information may further represent at least a section 
of a target and a length of at least one primer or probe complimentary to the target. 
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The hybridization thermodynamics may be calculated for a plurality 
of primers or probes complimentary to the target. 

The hybridization information may represents at least a section of a 
target and a primer or probe. 

5 A length of the target may be longer than a length of the primer or 

probe and wherein the hybridization thermodynamics are calculated for a best 
target/primer or target/probe complex and for competitive mismatch complexes. 

Hybridization information may represent at least a section of a target 
and a primer or probe and wherein a length of a target is longer than the length of 
10 the primer or probe and wherein the hybridization thermodynamics are calculated 
for a best target/primer or target/probe complex and for competitive target/primer 
or target/probe complexes. 

The method may further include calculating concentration of each 
species in a solution at a plurality of temperatures. 

15 Hybridization information may also represent a primer or probe and 

wherein the length of the target is longer than a length of the primer or probe and 
wherein the hybridization thermodynamics are calculated for a best target/primer 
or target/probe complex and for competitive mismatch complexes and wherein the 
method may further comprise calculating concentration of every species in a 

20 solution at a plurality of temperatures. 

The hybridization thermodynamics may be calculated for at least two 
best target/primer or target/probe complexes and for their corresponding 
competitive mismatch complexes and wherein the method may further comprise 
correcting for any interactions between the at least two best target/primer or 
25 target/probe complexes and tiieir components. 
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Further in carrying out the above objects and other objects of the 
present invention, a system for predicting nucleic acid hybridization 
thermodynamics is provided. The system includes a database of thermodynamics 
parameters, means for receiving hybridization information which represents at least 
5 one sequence, and means for receiving correction data. The system further includes 
receiving a first set of data which represents hybridization conditions, and means 
for calculating hybridization thermodynamics including net hybridization 
thermodynamics based on the hybridization information, the thermodynamic 
parameters, the correction data and the first set of data. 

10 The hybridization thermodynamics of individual single stranded, 

bimolecular and higher order complexes may be statistically weighted in a numerical 
process and the equilibrium concentration of each species is output. 

The correction data may include folding correction data and/or linear 
correction data. 

15 The thermodynamic parameters may include DNA thermodynamic 

parameters such as dangling end parameters. 

The DNA thermodynamic parameters may include coaxial stacking 
parameters and/or terminal mismatch parameters. 

The thermodynamic parameters may include RNA thermodynamic 
20 parameters and/or hybrid DNA/RNA thermodynamic parameters. 

The thermodynamic parameters may further include DNA loop 
thermodynamic parameters. 

The hybridization information may represent top and bottom strand 
sequences which form a duplex and wherein the hybridization thermodynamics are 
25 calculated for the duplex. 
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The hybridization information may also represent at least a section 
of a target and a length of at least one primer or probe complimentary to the target. 

The hybridization thermodynamics may be calculated for a plm-ality 
of primers or probes complimentary to the target. 

5 The hybridization information may represent at least a section of a 

target and a primer or probe. 

A length of the target may be longer than a length of the primer or 
probe and wherein the hybridization thermodynamics are calculated for a best 
target/primer or target/probe complex and for competitive mismatch complexes. 

10 Hybridization information may represent at least a section of a target 

and a primer or probe and wherein a length of a target is longer than the length of 
the primer or probe and wherein the hybridization thermodynamics are calculated 
for a best target/primer or target/probe complex and for competitive target/primer 
or target/probe complexes. 

15 The system may further include means for calculating concentration 

of each species in a solution at a plurality of temperatures. 

Hybridization information may also represent a primer or probe and 
wherein the length of the target is longer than a length of the primer or probe and 
wherein the hybridization thermodynamics are calculated for a best target/primer 
20 or target/probe complex and for competitive mismatch complexes and wherein the 
system may further comprise means for calculating concentration of every species 
in a solution at a plurality of temperatures. 

The hybridization thermodynamics may be calculated for at least two 
best target/primer or target/probe complexes and for their corresponding 
25 competitive mismatch complexes and wherein the system may further comprise 
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means for correcting for any interactions between the at least two best target/primer 
or target/probe complexes and their components. 

Still further in carrying out the above objects and other objects of the 
present invention, a computer-readable storage medium having stored therein a 
5 database of thermodynamics parameters and a computer program are provided. The 
computer program executes the steps of: a) receiving hybridization information 
which represents at least one sequence; b) receiving correction data; c) receiving a 
first set of data which represents hybridization conditions; and d) calculating 
hybridization thermodynamics based including net hybridization thermodynamics 
10 based on the hybridization information, the thermodynamic parameters, the 
correction data and the first set of data. 

The hybridization thermodynamics of individual single stranded, 
bimolecular and higher order complexes may be statistically weighted in a numerical 
process and the equilibrium concentration of each species is output. 

15 . The correction data may include folding correction data and/or linear 

correction data. 

The thermodynamic parameters may include DNA thermodynamic 

parameters. 

The DNA thermodynamic parameters may include dangling end 
20 parameters and/or coaxial stacking parameters. 

The DNA thermodynamic parameters may further include terminal 
mismatch parameters. 

The thermodynamic parameters may include RNA thermodynamic 
parameters and/or hybrid DNA/RNA thermodynamic parameters. 
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The thermodynamic parameters may further include DNA loop 
thermodynamic parameters. 

The hybridization information may represent top and bottom strand 
sequences which form a duplex and wherein the hybridization thermodynamics are 
5 calculated for the duplex. 

The hybridization information may represent at least a section of a 
target and a length of at least one primer or probe complimentary to the target. 

The hybridization thermodynamics may be calculated for a plurality 
of primers or probes complimentary to the target. 

10 The hybridization information may represent at least a section of a 

target and a primer or probe. 

A length of the target may be longer than a length of the primer or 
. probe and wherein the hybridization thermodynamics are calculated for a best 
target/primer or target/probe complex and for competitive mismatch complexes. 

15 Hybridization information may represent at least a section of a target 

and a primer or probe and wherein a length of a target is longer than the length of 
the primer or probe and wherein the hybridization thermodynamics are calculated 
for a best target/primer or target/probe complex and for competitive target/primer 
or target/probe complexes. 

20 The program may further execute the step of calculating 

concentration of each species in a solution at a plurality of temperatures. 

Hybridization information may also represent a primer or probe and 
wherein the length of the target is longer than a length of the primer or probe and 
wherein the hybridization thermodynamics are calculated for a best target/primer 
25 or target/probe complex and for competitive mismatch complexes and wherein the 

AO- 
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program may execute the step of calculating concentration of every species in a 
solution at a plurality of temperatures. 

The hybridization thermodynamics may be calculated for at least two 
best target/primer or target/probe complexes and for their corresponding 
5 competitive mismatch complexes and wherein the program may execute the step of 
correcting for any interactions between the at least two best target/primer or 
target/probe complexes and their components. 

The above objects and other objects, features, and advantages of the 
present invention are readily apparent from the following detailed description of the 
10 best mode for carrying out the invention when taken in connection with the 
accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIGURE 1 is a schematic drawing wherein multiple equilibria are 
considered for concentration calculations; 

15 FIGURE 2a is a schematic drawing of a user input interface wherein 

the user provides various input iaformation for a fibrst module of the invention; 

FIGURE 2b is a schematic drawing of a user output interface wherein 
a computer provides output information corresponding to the input information of 
Figure 2a; 

20 FIGURE 3a is a schematic drawing of a user input interface wherein 

the user provides various input information for a second module of the invention; 

FIGURE 3b is a schematic drawing of a user output interface wherein 
a computer provides output information corresponding to the input information of 
Figure 3a; 
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FIGURE 4a is a schematic drawing of a user input interface wherein 
the user provides various input information for a third module of the invention; 

FIGURE 4b is a schematic drawing of a user output interface wherein 
a computer provides output information corresponding to the input information of 
5 Figure 4a; 

FIGURE 5a is a schematic drawing of a user input interface wherein 
the user provides various input information for a fifth module of the invention; 

FIGURE 5b is a schematic drawing of a user output interface wherein 
a computer provides output information corresponding to the input information of 
10 Figure 5a; 

FIGURE 6 is a block diagram flow chart illustrating the solution of 
conservation equations of the present invention; 

FIGURE 7 is a schematic diagram illustrating multiplex PGR design; 

FIGURE 8 shows prediction of molecular beacon net hybridization 
15 thermodynamics; 

FIGURE 9 shows simulation of molecular beacon hybridization 
concentrations at temperatures from 0 to 100 °C; 

FIGURE 10 is a diagram of match vs. mismatch hybridization; 

FIGURE 11 shows match vs. mismatch hybridization simulation at 
20 different temperatures; 

FIGURE 12 shows a general case of competitive hybridization 
equilibria that can be solved using the described numerical methods; and 
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FIGURE 13 is an example of simultaneous equations for the general 
five molecule case. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

In general, the method and system of the present invention include 
5 rigorous thermodynamic parameterization for Watson-Crick base pairs, internal 
mismatches, terminal mismatches, terminal dangling ends, co-axial stacking 
interactions, sodium and magnesium salt dependence, denaturants (urea, formamide, 
DMSO). In addition, loop parameters for hairpins, internal loops, bulges, and 
multibranched loops are included. For DNA essentially all the parameters have 

10 been previously published or all included in the Appendix hereto. Specifically, the 
parameters which have been published include Watson-Crick parameters, sodium 
dependence, GT, GA, CT, AC, AA, CC, GG, and TT mismatches. The 
parameters included hereiu include dangliug ends, terminal mismatches, DNA loop 
parameters, and co-axial stacking parameters. For RNA, the parameters have been 

15 published by Douglas H, Turner et al. For DNA/RNA hybrid duplexes, the 
parameters have been published by Naoki Sugimoto. 

The method and system are adapted for fiiture implementation of 
parameters for modified nucleosides (including but not limited to inosine, 5- 
nitroindole, PNA, MOE-modified RNA, and iso-bases). With these parameters, it 

20 is possible to predict the melting temperature, Tm, of a duplex within 2°C on 
average. Correction for surface effects for DNA chip arrays is also implemented. 
In addition to predicting duplex hybridization, the software accounts for single- 
strand secondary structure. This is accomplished by a new numerical procedure for 
solving complex coupled equilibria (multi-state model). With this approach, it is 

25 possible to accurately predict not only the Tm for hybridization but also the 
concentration of every species in the solution (e.g. match duplex, mismatch, duplex, 
folded target, folded primer, primer dimer, etc.) at every temperature from 0 to 
100°C. Thus, it is possible to use this software to design oligonucleotide 
hybridization with optimized temperature, salt, and strand concentrations. 
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PredictinQ Accurately Primer Target Interaction Stability 

The stability of a primer/target or probe/target complex can be 
described by the free energy of association of the probe and the target. The most 
accurate way to calculate free energy of association is to use the nearest-neighbor 
5 model with accurate thermodynamic parameters. Thermodynamic parameters 
should account for Watson-Crick base pairs (SantaLucia et al., 1996; AUawi & 
SantaLucia, 1997; SantaLucia, 1998), single mismatches (AUawi & SantaLucia, 
1997), terminal mismatches (disclosed herein), dangling ends (Bommarito, Pugret 
& SantaLucia, 2000) and possibly double mismatches. Proper calculation of the 

10 monovalent and divalent salt dependence is also important (SantaLucia, 1998). 
Other loop motifs for hairpins, bulge, internal loops and multi-branched loops are 
important for single strand secondary structure prediction, but are often very 
crudely approximated. Moreover, when primer and target folding can occur, a set 
of coupled equilibria should be used to model the system. The nearest-neighbor 

15 model needs to be used to determine the equilibrium constant of each equilibrium. 
The determination of possible primer or target folding can be addressed by using 
secondary-structure prediction algorithms like M. Zuker's MFOLD (Zuker, 1989). 

Secondary Structure and Net Hybridization Thermodynamics 

Species Concentration Calculations 

20 Consider a system of strands SI and 82 with four states: folded 

target, folded probe, probe bound to target, and random coil target and probe. The 
model can be described by three equilibria as shown in Figure 1. 

The concentrations of every species for such a system can be 
analytically determined. The three equilibrium constants for such a system are 
25 shown below: 

S1+S2=DH = (I) 

-14- 
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Sl^Hl = ^ 



S2^H2 - ^ (3) 



where SI, S2, HI, H2, and DH are the random coil SI, the random coil S2, the 
folded strand HI, the folded strand H2, and the double helix DH, respectively. The 
5 conservation of SI and S2 leads to the following equations: 

C^°"" = Sl+ m+ DH (4) 

Cjr' = S2+H2+ DH (5) 

where d^?^' are the total concentrations of SI and S2. [DH] and [S2] can be 
10 expressed as a function of [SI] by substituting the [HI] obtained from Equation 2, 
in Equation 6, and then substituting the [DH] obtained by Equation 6 in Equation 
1. 

[DH] = Cjf"^ - [SI] - K^[S1] (6) 



^S2]= ^s2]^C^r'-[Sl]-K,[Sl] 



15 Substitution of [H2], [DH], and [82] from Equations 3, 6 and 7 in 

Equation 5 leads to an expression of Ki that can be rearranged as a quadratic 
equation in [SI]: 

[Slf (K, + K,) + [S1](K, Cs2 + K3 + K2 K3 - K, Csi + + 1) - (K3 + 1) Cs, = 0 

(8) 

20 This equation is simplified by making the following substitutions: 
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a = (Ki + K2 K,) (9) 
b = (K, C^r' + K3 + K2K3-K1 Clr' + K2 + 1) (10) 
c = (K3 + 1) (11) 

The physical solution of the quadratic equation (i.e. positive root) is (Press, 1999): 

ISU. -'*^^ (12) 



or 



^ - b- -Jb^ - 4ac 



The second equation has better numerical stability (Press, 1999). 
[DH], [S2], [HI], and [H2] can then be calculated using Equations 1-3. 
10 Determination of Net Free Ener^ 

The net free energy of hybridization is calculated as follows: 



^G\,= -RT]nK, (14) 



where 



^ [DH] ^^^^ 

\.^^ single stranded ^i^^ single stranded^ 



15 where [SI single stranded] and [S2 single stranded] are the concentrations of SI and 
S2 either in the random coil state or the hairpin states, at the temperature of the 
simulation. Using the conservation of SI and S2, Equation 14 is rewritten as 
follows: 
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AG%,, ^^RT In ^ [Dm)^lT' ' W^]) ^^^^ 

Note that AG^j^ has the unusual property that it depends on the total strand 
concentrations, C^°^^ and (Jsi^ . The net free energy expresses the duplex formation 
equilibrium free energy corrected for secondary-structure formation in the single 
5 strands. 



Determination of Net Melting Temperature 

If the strands are non self-complementary two cases have to be 
considered depending on the relative strand concentrations: 

1) If SI is the limiting reagent (C^^ < ^^ii"). at Tm: 



10 {DH]=~Clr' (17) 

The concentrations of strands [SI] and [S2] are given by the following relations: 



= 2^sT' + + [^S'l] (18a) 



^ Total 



Cjr' = 5Cjr' + J£3[S2]+[S2] (19a) 



^ Total _ X./^^otal 

15 [^2]= " ~^ (19b) 

+ 1 
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The replacement of [SI] and [S2] in Equation 1 gives: 

' (20a) 



/nr Total J_ /nr Tbto/ 
^52 2 



0 = 1(^7''^"^/ _ + + -^2^3 + -^3 + 1 (20b) 

2 ATi 



Using the relation AG°t= -R T In K, Equation 20 is arranged as follows: 



1 f -A(7°r(2)-AGV(3) -AGV(3) A AGV(l) 



2 



(21) 

AG°t can then be decomposed as AG°x=AH*'- T AS° (assuming AC p^o) to obtain: 

-A//^(2)-A//°(3)+A//°(l) A5'*(2)-t-A5"(3)-A5°(l) -A/Z'^CSHA/Z^Cl) A5°(3)-A^°(l) A//°(l) -A5°(l) 

(22) 

10 The above equation can be solved by bisection or other numerical techniques to find 
T. This solution is the net melting temperature. 

2) If S2 is the limiting reagent (ClT^ < CFsT\ the following relation can 
be deduced by a similar approach: 



15 Again, application of the bisection method to an equation symmetric to Equation 22 
affords the net melting temperature. 
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If the strand S is self complementary, the reactions are described by 
the following equilibria: 



S + S -DH = ^ (24) 



S^H K,= ^-^ (25) 



AtT^iiDHl-^Cl""" (26) 
The strand conservation equation is: 

Cf = [5] + [/Tl + 2[DH] (27) 
Insertion of [H] and [DH] from Equations 25 and 26 in Equation 27 leads to: 



Cj*-' = 2{K^ + 1)[5] o [5-] = (28) 



10 Introduction of [S] in Equation 24 gives: 



= ^= 0= + 2i., - K,cr . 1 (29) 



Using the relation AG°t= -R T In K, Equation 29 is rearranged as follows: 



-2AG°r(2) -AGV(2) AGV(l) 

0=e +2e -Cf^'e ''^ +1 (30) 
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AG°T can then be decomposed as AG''t=AH°- T AS° to obtain: 

-2A/f°(2) 4-2A>S'°(2) -Ai/°(2) 'fA5°(2) -Aj^°(l) +A^°(1) 

0=e e ^ +2e e ^ - Cf '^'e ^ ^ +1 (31) 

This equation can be solved by bisection to afford the net melting temperature. 

An experimentally validated example of the accuracy of the net 
5 hybridization thermodynamics is shown in Figure 8 for molecular beacons. At the 
top of Figure 8 are the predicted thermodynamics for simple duplex formation 
assuming no competing single strand secondary structure. Using Module 1 of the 
invention, these results are similar to what would be predicted using other 
commercial software (such as oligo 6.0), though our thermodynamic database 

10 includes the dangling end effects and salt corrections are more accurate than other 
software. The middle of Figure 8 shows the single strand folding at the molecular 
beacon as output firom DNA-MFOLD. The bottom table of Figure 8 shows the 
experimentally determined A6 (effective) and Tm (effective) published in Bonnet et 
al. 1999, as well as the effective Tm and A6 (effective) predicted with Module 1 

15 using the coupled equilibria calculations. Note the close agreement between 
experiments and predictions in the bottom table and the disagreement between 
experiments and the predictions using the naive simple hybridization calculation (top 
table of Figure 8). Also note the good agreement in the bottom table for the fully 
matched A-T sequence and mismatch A-A, A-C, and A-6 sequences, thus validating 

20 the mismatch parameters. 

Further, the net hybridization calculations can be extended to 
different temperatures as shown in Figure 9, to reveal how the concentrations of all 
species change with temperature. Given the extinction coefficients and fluorescence 
quantum yields, the concentration vs. temperature profiles shown in Figure 9 can 
25 be used to calculate the fluorescence vs. temperature profile (not shown), thereby 
allowing the prediction of the temperature which produces the maximum 
fluorescence signal and minimum background fluorescence signal. 
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Another manifestation of the concentration calculations is for match 
vs. mismatch discrimination (Figure 10), whereby the concentrations of all species 
at all temperatures can be calculated (Figure 11). For the particular case shown, 
optimal match vs, mismatch discrimination is predicted to occur at C^C. The 
5 concentration calculations can be generalized for cases in which molecules can form 
many different competing unitnolecular, biomolecular, and higher order complexes 
(Figure 12) using generalized equations such as shown in Figure 13 for the five 
molecule case, and solved using the algorithm in Figure 7. 

Alsorithm 

10 The hybridization prediction algorithm of the present invention is 

based on a nearest-neighbor-model analysis of the sequences. The algorithm 
accounts for structural motifs including Watson-Crick base pairs (AUawi and 
SantaLucia, 1997; SantaLucia, 1998; Sugimoto et al., 1995; Xia et al., 1998), 
single internal mismatches (AUawi and SantaLucia, 1997; Allawi and SantaLucia, 

15 1998; Allawi and SantaLucia, 1998; Allawi and SantaLucia, 1998; Kierzek et al., 
1999; Peyret et al., 1999; SantaLucia, 1998), double mismatches (Allawi and 
SantaLucia, 1997) coaxial-stacking interfaces (disclosed herein) (Walter and Turner, 
1994), terminal mismatches (disclosed herein) (Freier et al., 1986) and dangling 
ends (Bommarito et al., 2000; Freier et al., 1986). Once the motifs are identified 

20 and their thermodynamic contributions are added, the sum may be corrected for salt 
effects (sodium and magnesium) and the net hybridization is calculated when 
appropriate. 

Algorithm Functions 

A first or main module of the algorithm calculates the hybridization 
25 thermodynamics (AH°, AS°, AG^37, T^) of a given duplex. Net hybridization 
accounting for secondary structure in both strands is also calculated. 
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Parameterization 

Parameters are organized in three arrays. The first array contains 
internal element parameters: Watson-Crick nearest neighbors and single mismatch 
nearest neighbors. The second array contains terminal element parameters: terminal 
5 mismatches and dangling ends. A single parameter is used to account for double 
mismatches except for tandem G T mismatches, which are explicitly enumerated 
(AUawi & SantaLucia, 1997). The third array contains coaxial-stacking parameters 
(contained herein). 

For DNA sequences, the thermodynamic contribution of all Watson- 
10 Crick nearest neighbors and single internal mismatches has been systematically 
studied (AUawi and SantaLucia, 1997; Allawi and SantaLucia, 1998; AUawi and 
SantaLucia, 1998; Allawi and SantaLucia, 1998; Peyret et al., 1999). A limited 
number of sequences containing double mismatches has also been studied (Allawi 
and SantaLucia, 1997). The contributions of dangling ends (Bommarito et al., 
15 2000) have also been systematically analyzed. Salt corrections are available for 
sodium in the range 0.01 to 1 M (SantaLucia, 1998). 

For RNA sequences, the thermodynamic contribution of all Watson- 
Crick nearest neighbors has been systematically studied (Xia et al., 1998). A limited 
number of sequences contairiing single mismatches has also been studied (Kierzek 
20 et al., 1999). The contribution of dangling ends and terminal mismatches has also 
been systematically analyzed (Freier et al., 1986). No salt correction has been 
developed for RNA and therefore the DNA salt corrections are assumed. These 
corrections are likely to be deficient in the case of RNA. 

For DNA/RNA hybrids, the thermodynamic contribution of all 
25 Watson-Crick nearest neighbors has been systematically studied as well as a limited 
number of sequences containing single mismatches (Sugimoto et al., 1995). As no 
salt correction has been developed for DNA/RNA hybrids, the DNA corrections are 
assumed. The applicabUity of these corrections to DNA/RNA hybrids has not been 
tested. The parameter arrays are designed to easily accommodate implementation 
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of new parameters and salt corrections including thermodynamics parameters for 
modified bases and denaturant effects. 

Correction for Hybridization to DNA Microchips 

A linear correction of the free energy is implemented in the algorithm 
5 of the invention to correct for hybridization to DNA microchips: 

AG°37(microchip) = aAG°37(solution) +b (32) 

where a and b are user defined real coefficients. Fotin et al. (Fotin et al., 1998) 
showed that a linear relationship could be used to relate the free energies obtained 
for hybridization in solution and on microchip surfaces. However, the relation 
10 between thermodynamics measured in solution and thermodynamics measured using 
microarrays is still unclear and appears to be different depending on the 
manufacture and type of microarrays. 

User Interface: Input and Output 

Figure 2a shows the user interface input. The users enter the 
15 sequence of each strand, the hybridization conditions (hybridization temperature, 
strand concentrations, and monovalent cations and concentrations), and 
thermodynamic corrections for single strand folding. Figure 2b shows the output 
corresponding to the input in Figure 2a. 

The algorithm can be used via the Internet at: 
20 http : ll\ si 1 . chem. wavne . edu/Hvther/hvtherm Imain. html . The algorithm may be 
written in FORTRAN 77 and run on UNIX environment or other languages and 
environments. 
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Molecular Beacons 

The algorithm may be used to predict the thermodynamics of a set 
of literature measurements for molecular beacons (Bonnet et al., 1999). Molecular 
beacons are high specificity probes that are efficient for mutation analysis 
5 (Giensendorf et al., 1998) and multiplex detection of single nucleotide variations 
(Marras et al., 1999). The design and efficiency optimization of these beacons is 
helped by hybridization thermodynamics prediction. Bonnet et al. studied, the 
hybridization of the molecular beacon ^ CGC, TCC, CAA, AAA, AAA, AAA, 
CCG AGC G^' to a set of four different targets including a perfect match duplex, 
10 and three different duplexes containing one mismatch. Free energy and enthalpy for 
duplex folding may be calculated using the DNA MFOLD program 
(http://mfold2.wustl.edu/ — mfold/dna/forml.cgi). These parameters may then 
incorporated as secondary structure corrections in Figure 2a. 

The software to implement the algorithm may be written in 
15 FORTRAN, C-^"^, Visual Basic, HTML, and JAVA script computer languages. 
Two graphical user interfaces may be provided: Windows application and web 
browser format. The software may run on IBM/PC, Sun, and Silicon Graphics 
platforms. 

The software may be written in several modules as described below. 

20 A. Interactive Mode: Command Line Interface in MS-DOS 

MODULE 1 (As Previouslv Described above) 

Function . Module 1 predicts the hybridization thermodynamics of a given 
duplex (DNA/DNA, RNA/RNA, or DNA/RNA). 
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Input (Figure 2a) 

Input of Sequences 

1 . Only the following characters are accepted: A, a, C, c, G, g, 
T, t, U, u, /,*,+. Single blank characters and numbers will be automatically 

5 edited, but more than one carriage return is not permitted. 

2. If the duplex contains a dangling end on a strand, the 
sequence of the other strand should contain a * at the corresponding position. (This 
is very important to include for primer binding to a large target sequence). Note: 
The top strand must be entered in 5' to 3' orientation, but the bottom strand must 

10 be entered in 3' to 5' orientation. Also, a " + " must be added at the end of each 
sequence. There is a length limit of 1024 characters for sequence entries. In 
module 1, it is important to be sure that both sequences have the same length. 

Example: AAAACCCCTGA + 
*TTTGGGGAC* + 

15 3. Only the bottom strand may contain coaxially stacked 

nucleotides. A should be inserted at the site of a strand nick {i.e. between the 
coaxially stacked nucleotides). This feature is useful for predicting stacked 
hybridization stability. 

Example: AAAACCCCC + 
20 TTTT/GGGG+ 

Input of Salt and Strand Concentrations 

The monovalent salt should be the sum of all monovalent cation 
concentrations in a solution in units of molarity. For example, a solution of 100 
mM KCl, 50 mM NaCl, 10 mM Na2P04, 0. 1 mM Na2EDTA would account for a 
25 total of 0. 1702 M monovalent. The thermodynamic predictions are applicable over 
a salt range of 0.01 to 1 M monovalent cation. The correction applied is from 
SantaLucia (1998) Proc. Natl Acad. Sci. 95, 1460. The sodium correction applies 
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for oligonucleotides with fewer than about 30 base pairs. For longer duplexes a 
polymer correction is required, but this is not currently implenoiented. 

Strand concentrations are entered in units of molarity. The program 
will accept virtually any physically relevant strand concentration. 

5 Hybridization temperature is in Celsius degrees. The limits are 0 to 

100 degrees. 

Special corrections for single-stranded secondary structure and for 
surface corrections for hybridization arrays can be input. The units for input AG° 
are kcal/mol. To determine estimates of single-strand folding energies, see Michael 
10 Zuker's RNA or DNA-MFOLD servers (see http: 
//mfold2.wsutl.edu/-'mfold/dna/formLcgi). The current thermodynamic prediction 
software incorporates the special corrections for single-stranded secondary structure 
and for surface corrections for hybridization arrays. 

For DNA chip arrays, a linear correction can be applied. The user 
15 inputs the slope and intercept coefficients. Based on the work of Mirzabekov 
group, a slope of +1.1 and intercept of +3.2 are appropriate (see Fotin et al. 
(1998) Nucleic Acids Res. 26, 1515-1521). 

Output (Figure 2b) 

Module 1 outputs the hybridization thermodynamics at 1.0 M NaCl 
20 and 37 ""C (the conditions under which the thermodynamic predictions are most 
accurate), under the salt temperature conditions specified by the user, and also 
displays the net hybridization Tm and AG° if the user specifies that special 
corrections are needed (this allows for single-strand secondary structure of both the 
target and probe DNA to be accounted and for surface effects of chip arrays). 
25 Predictions of AG° , AH ° , AS ° , and Tm are provided. 



-26- 



wo 01/94611 



PCT/USOl/18424 



MODULE 2 

Function . Module 2 finds the best primers of given length 
complementary to a long target nucleic acid. DNA/DNA, RNA/RNA, DNA/RNA 
hybridization types are accepted. The user selects the number of primers to output, 
5 and the program finds the most stable primers and gives their hybridization position 
and thermodynamics of each primer. 

Input (Figure 3a) 

Input of Salt and Strand Concentrations 

The input of strand and salt concentrations is similar to Module 1. 

10 Input of Sequences 

The target sequence is input as in Module 1. 

Output (Fisure 3b) 

Primer Length and Number of Best Primers 

Module 2 displays "number of best primers" best primers of length 
15 "primer length" in order of decreasing stability. 

Output 

Module 2 outputs "number of best primers" best primers of length 
"primer length" in order of decreasing stability along with their hybridization 
thermodynamics . 

20 MODULE 3 

Function , Module 3 walks a given primer along a given target and 
finds the thermodynamics for the best target/primer complex and for the competitive 
target/primer complexes: DNA/DNA, RNA/RNA, DNA/RNA, hybridization types 
are accepted. 
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Input (Figure 4a) 

Input of Sequences 

The input is similar to Module 1. The target has to be longer than 

the primer. 

5 Input of Salt and Strand Concentrations 

The input of salt and strand concentrations is similar to Module 1 . 

Percent Stability p of Alternative Binding Sites 
Compared to the Most Stable Bindmg Site 

This parameter excludes all competitive sites that are not within the 
10 defined percent of the best primer stability. If the best primer stability is -5 
kcal/mol and p = 10 then any competitive site of energy higher than -5 + (10/100*5) 
= -4.5 kcal/mol will not be displayed. 

Number of Base Pairs Required to Compute the Solution 
This parameter excludes all competitive sites that contain less 
15 Watson-Crick base pairs than the defined value. 

Output (Figure 4b) 

Module 3 outputs the best primer binding site and the competitive 
binding sites that pass the filtering criteria (percent stability p of alternative binding 
sites compared to the most stable binding site and number of best primers). 

20 MODULE 4 

Function . Batch mode calculations (see below). 
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MODULE 5 

Function . Module 5 is a combination of Modules 2 and 3 and finds 
the n best primers of given length complementary to a given section of a target and 
display the thermodynamics of the target/primer system(s). Then, each best primer 
5 is walked along the whole target to find the competitive hybridization sites. The 
thermodynamics of the target/primer systems at these alternative sites is then 
displayed. DNA/DNA, RNA/RNA, DNA/RNA, hybridization types are accepted. 

Input (Figure 5a) 

Input of Sequences 
10 The target sequence is input as in Module 1. 

Input of Salt and Strand Concentrations 

The input of salt and strand concentrations is similar to Module 1 . 

Sequence Section Where to Find the Best Primers 
Module 5 finds the best primers in the target region ranking from 
15 "position of initial nucleotide" to "position of final nucleotide" . Note that Module 
5 then looks for competitive sites of each best primers in the whole target. 

Percent Stability of Alternative Binding Sites 
Compared to the Most Stable Binding Site 

The function of this parameter is the same as in Module 3. This 
20 parameter is input for each best primer corresponding to the "number of best 
primer" specified. 

Number of Base Pairs Required to Compute the Solution 
The function of this parameter is the same as in Module 3. This 
parameter is input for each best prhner corresponding to the "number of best 
25 primer" specified. 
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Output (Fisure 5b) 

Primer Length and Number of Best Primers 

Module 5 displays "number of best primers" best primers of length 
"primer length" by order of decreasing stability. 

5 Output 

Module 5 displays "number of best primers" best primers and their 
competitive sites by order of stability along with their hybridization 
thermodynamics. The best primer and its ranked competitive hybridization sites are 
listed first. Then, the second best primer is listed with its competitive hybridization 
10 sites. 



MODULE 6 



Function . Module 6 is similar to Module 3 and walks a given primer 
along a given target and finds the thermodynamics for the best target/primer 
complex and for the competitive target/primer complexes: DNA/DNA, RNA/RNA, 
15 DNA/RNA, hybridization types are accepted. Then, Module 6 simulates the 
concentration of every species at every degree from 1 to 100 °C, as illustrated in 
Figure 6. 

Input (Not Shown) 

Input of Sequences 

20 The input is sunilar to Module 1. The target has to be longer than 

the primer. 

Input of Salt and Strand Concentrations 

The input of salt and strand concentrations is similar to Module 1 . 
Percent Stability p of Alternative Binding Sites 
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Compared to the Most Stable Binding Site 

This parameter excludes all competitive sites that are not within the 
defined percent of the best primer stability. If the best primer stability is -5 
kcal/mol and p = 10, then any competitive site of energy higher than: -5 + 
5 (10/100*5) = -4,5 kcal/mol will not be displayed. 

Number of Base Pairs Required to Compute the Solution 
This parameter excludes all competitive sites that contain less 
Watson-Crick base pairs than the defined value. 

Correction for Target/Target Interaction, Target folding, 
10 Primer/Primer Interaction and Primer Folding 

The user is asked if he wants to correct for the interactions above. 

If the answer is "y", the user is prompted for AH°37 corresponding to the 

interaction. Secondary structure thermodynamics can be determined using the 

Zuker algorithm as discussed in Module 1 section. 

15 Output (Not Shown) 

Concentration Output Filename 

The results from the concentration simulations (concentration of 
species at every temperature) are saved in this file. 

Output 

20 Module 6 outputs the best primer binding site and the competitive 

binding sites that pass the filtering criteria (percent stability p of alternative binding 
sites compared to the most stable binding site and number of best primers). The 
concentration simulations are saved in a file specified by the user. 
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MODULE 7 

Function , Module 7 is a combination of Modules 2 and 5 and finds 
the n best primers of given length complementary to a given section of a target and 
display the thermodynamics of the target/primer system(s). Then, each best primer 
5 is walked along the whole target to find the competitive hybridization sites. The 
thermodynamics of the target/primer systems at these alternative sites is then 
displayed. DNA/DNA, RNA/RNA, DNA/RNA hybridization types are accepted. 
Then, Module 7, like Module 6, simulates the concentration of every species at 
every degree from 1 to 100°C, as illustrated in Figure 6. 

10 Input (Not Shown) 

Input of Sequences 

The target sequence is input as in Module 1 . 
Input of Salt and Strand Concentrations 

The input of salt and strand concentrations is similar to Module 1 . 

15 Sequence Section Where to Find the Best Primers 

Module 7 finds best primers in the target region ranking from 
"position of initial nucleotide" to "position of final nucleotide. " Note that Module 
7 then looks for competitive sites of each best primers in the whole target. 

Percent Stability of Alternative Binding Sites 
20 Compared to the Most Stable Binding Site 

The function of this parameter is the same as in Module 3. This 

parameter is input for each best primer corresponding to the "number of best 

primer" specified. 
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Number of Base Pairs Required to Compute the Solution 

The function of this parameter is the same as in Module 3. This 

parameter is input for each best primer corresponding to the "number of best 

primer" specified. 

5 Correction for Target/Target Interaction, Target folding, 
Primer/Primer Interaction and Primer Folding 

For each best primer, the user is asked if he wants to correct for the 

interactions above. If the answer is **y ", the user is prompted for AH° and AG°37 

corresponding to the interaction. Secondary structure thermodynamics can be 

10 determined using the Zuker algorithm as discussed in Module 1 section. 

Concentration Output Filenames 

For each best primer, the results from the concentration simulations 
(concentration of species at every temperature) are saved in this file. The user has 
to select a different filename for each best primer. 

15 Output (Not Shown) 

Output 

Primer Length and Number of Best Prhners 

Module 7 displays "number of best primers" best primers of length 
"primer length" by order of decreasing stability. 

20 Module 7 displays "number of best primers" best primers and their 

competitive sites by order of stability along with their hybridization 
thermodynamics. The best primer and its ranked competitive hybridization sites are 
listed first. Then, the second best primer is listed with its competitive hybridization 
sites. For each best primer, a file named by the user contains the concentration 

25 simulations. 
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Module 7 allows Hie user to design optimal primers for applications 
where multiple simultaneous hybridization reactions are occurring, including match 
vs. mismatch hybridization, molecular beacons, DNA oligonucleotide arrays, and 
multiplex PGR. 

5 One commercially important example for the use of Module 7 for 

primer design in a complex hybridization solution is Multiplex PGR, as shown in 
Figure 7. Module 7 allows the user to design optimal primers for Multiplex PGR 
where multiple primers have equal stabilities ia binding to the target DNA, Several 
primers must be designed to specifically bind to different sites on target DNA at a 
10 given temperature with minimal background binding to mismatch sites and with 
minimal cross-hybridization between pairs of primers. 

Module 7 rninhmzes potential primer dimer formation and mismatch 
hybridization for all combinations of input primers. Module 7 optimizes primer 
sequence position, length, and concentration for each primer in relation to all other 
15 species in solution and provides a hybridization profile at all temperatures from 0 
to 100°G. 

Batch Mode 

MODULE 4 

Function . Module 4 allows any of the previous modules to be run 
20 in batch mode using text files to submit the input and having the data output as text 
files also. 

Tvpe of Input Files 

There are two types of input files: 1) parameter input file, and 2) 
sequence input file. Parameter input files describe what modules to run with what 
25 hybridization parameters and on how many sequences to run them. Example of 
parameter input files for each module with comments are given in the "Batch mode 
parameter files folder. " Sequence files contain the sequences that are going to be 
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hybridized in the conditions described by the parameter input files. Examples of 
parameter input files for each module with comments are given in the "Batch mode 
sequence files folder." 

Note that a parameter file can successively run different modules on 
5 various different sequences. 

The user is successively asked for the names of the parameter input 
file, the sequence input file and the thermodynamic data output file. Note that these 
files have to be in the directory containing the executable version of the software. 
Ou^ut files will also be created in this same directory. Names of the concentration 
10 simulation files are specified in the parameter input files. 

Examples of Batch Mode Parameter Files 

Comments in parentheses describe the meaning of each entry (note 
that an actual parameter file must not contain these comments). 

DUP (Module 1: Simple duplex calculations) 

15 1 (Number of sequences to apply this parameter file to) 

1 (Monovalent cations concentration mol/L)) 

1 Mg^"^ concentration mol/L) 

37.0 (Hybridization temperature) 

4e-4 (Top strand concentration mol/L) 

20 4e-4 (Bottom strand concentration mol/L) 

1 (Correction for microchips: slope) 

0 (Correction for microchips: intercept) 

0 (Correction for top strand folding: AG°37) 

0 (Correction for top strand folding: AH°37) 

25 0 (Correction for bottom strand folding: AG^sy) 

0 (Correction for bottom strand folding: AH'^sy) 

END (End of file required) 
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NBP (Module 2: N-best primers) 

1 (Number of sequences to apply this parameter file to) 

1 (Monovalent cations concentration mol/L) 

1 (Mg^"^ concentration mol/L) 

5 37 (Hybridization temperature) 

4e-4 (Top strand concentration mol/L) 

4e-4 (Bottom strand concentration mol/L) 

4 (Primer Length) 

3 (Number of best primers) 

10 1 (Correction for microchips: slope) 

0 (Correction for microchips: intercept) 

END (End of file required) 



PWA (Module 3 primer walk match vs. mismatch sites identification) 

1 (Number of sequences to apply this parameter file to) 

15 1 (Monovalent cations concentration mol/L) 

1 (Mg^"^ concentration mol/L) 
37 (Hybridization temperature) 
4e-4 (Top strand concentration mol/L) 
4e-4 (Bottom strand concentration mol/L) 

20 90 (Percent window of best primer stability for alternative sites) 

2 (Number of WC base pairs required to compute the solution) 
cgcg+ (Primer sequence, + required) 

1 (Correction for microchips: slope) 

0 (Correction for microchips: intercept) 

25 END (End of file required) 



BPW (Module 5 displays "number of best primers" best primers and their 
competitive sites by order of stability along with their hybridization 
thermodynamics) 

1 (Number of sequences to apply this parameter file to) 
30 1 (Lower limit of primer search area) 

10 (Upper limit of primer search area) 
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1 Mg^"^ concentration mol/L) 

37 (Hybridization temperature) 

4e-4 (Top strand concentration mol/L) 

4e-4 (Bottom strand concentration mol/L) 

5 4 (Primer length) 

1 (Number of best primers) 

1 (Correction for microchips: slope) 

0 (Correction for microchips: intercept) 

800 (Percent window of best primer stability for alternative sites) 

10 2 (Number of WC base pairs required to compute the solution) 

END (End of file required) 

PWC (Module 6 primer walk with concentration calculations) 

1 (Number of sequences to apply this parameter file to) 
1 (Monovalent cations concentration mol/L) 

15 1 Mg^"*" concentration mol/L) 

37 (Hybridization temperature) 

4e~4 (Top strand concentration mol/L) 

4e-4 (Bottom strand concentration mol/L) 

90 (Percent window of best primer stability for alternative sites) 

20 2 (Number of WC base pairs required to compute the solution) 

cgcg+ (Primer sequence, + required) 

1 (Correction for microchips: slope) 

0 (Correction for microchips: intercept) 

0 (Correction for target folding: AC^gy) 

25 0 (Correction for target folding: AH°) 

0 (Correction for target/target interaction: AG°37) 

0 (Correction for target/target interaction: AH°) 

0 (Correction for primer folding: AG^jy) 

0 (Correction for primer folding: AH"") 

30 0 (Correction for primer/primer interaction: ACgy) 

0 (Correction for primer/primer interaction: AH"") 

outconc (Concentration ouq)ut file name) 
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END (End of file required) 





BWC 


(Module 7, N-best primers, primer walk, and concentration 






calculations) 




1 


(Number of sequences to apply this parameter file to) 


5 


1 


(Lower limit of primer search area) 




10 


(Upper limit of primer search area) 




1 


(Monovalent cations concentration mol/L) 




0 


Mg^^ concentration mol/L) 




37 


(Hybridization temperature) 


10 


4e-4 


(Top strand concentration mol/L) 




4e-4 


(Bottom strand concentration mol/L) 




4 


(Primer length) 




1 


(Number of best primers) 




1 


(Correction for microchips: slope) 


15 


0 


(Correction for microchips: intercept) 




800 


(Percent window of best primer stability for alternative sites) 




2 


(Number of WC base pairs required to compute the solution) 




0 


(Correction for target folding: AG^^-j) 




0 


(Correction for target folding: AH"") 


20 


0 


(Correction for target/target interaction: AG°37) 




0 


(Correction for target/target interaction: AH°) 




u 


^v^orrecuon ror pruner roiumg. avj 27) 




0 


(Correction for primer folding: AH^) 




0 


(Correction for primer/primer interaction: AG°37) 


25 


0 


(Correction for primer/primer interaction: AH°) 




outconc 


(Concentration output file name) 




END 


(End of file required) 




PPW 


(Module 8: walk a primer along itself to find interaction sites. PWA 






with probe = primer) 


30 


1 


(Number of sequences to apply this parameter file to) 




1 


(Monovalent cations concentration mol/L) 
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1 


Mg^^ concentration mol/L) 


37 


(Hybridization temperature) 


4e-4 


(Primer concentration mol/L) 


900 


(Percent window of best primer stability for alternative sites) 


2 


(Nmnber of WC base pairs required to compute the solution) 


END 


(End of file required) 



Examples of Batch Mode Sequence Files 

For Module 1 : Pup 
1 (Sequence number) 

10 agcgca+ (Top strand sequence) 

tcgcgt+ (Bottom strand sequence) 

For Module 2: NBP 

1 (Sequence number) 

agcgca+ (Target sequence) 

15 For Module 3: 

1 (Sequence number) 

cgcctgcggccc+ (Target sequence) 

For Module 5: bpw 
1 (Sequence nmnber) 

20 cgcctgcgccc+ (Target sequence) 

For Module 6: pwc 

1 (Sequence number) 

agcgca+ (Target sequence) 

For Module 7: bwc 
25 1 (Sequence number) 

agcgca-l- (Target sequence) 
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For Module 8: ppw 

1 (Sequence number) 

agcgca+ (Primer sequence) 



Example of Batch Mode Parameter and Sequertce Files 
5 to Run Different Modules Successively 



Parameter File : 

DUP (Executes Module 1) 

2 (Apply to Module 1 to 2 sequence sets) 

0.05 
10 1.5e-3 
37.0 
le-6 
2e-7 
1 

15 0. 

-2.12 
-37.3 
0 
0 

20 PWC (Executes Module 6) 

1 (Apply to Module 6 to 1 sequence set) 

0.16 

0.0025 

37 

25 lOe-9 
le-9 
800 
8 

TCGAACGTAC+ 
30 1 
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0 
0 
0 
0 

5 0 
0 
0 

0 
0 

10 outwash 

DUP (Executes Module 1) 



4 (Apply to Module 1 to 4 sequence sets) 

1 

0 

15 37.0 
le-6 
le-6 
1 
0 

20 0 
0 
0 
0 

END 

25 Other modules can be similarly appended. 
Sequence File 

1 (input for Module 1) 

ttgcctaggggaccaggtccaact + 
aacggatcccctggtccaggttga + 
30 2 (input for Module 1) 
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ttgcctaggggaccaggtccaact + 
aacggatcccctggtccaggttga + 
3 

C AGCTTGC ATGA AA AGCTTGCGTGT + 
5 4 

AAAAAA+ 
TTTTTT+ 
5 

acgcgc + 
10 tgcgcg+ 
6 

gggaaagggg + 
*cctttccc*+ 
7 

15 tttaaattt+ 
aaatttaaa+ 
8 

cgcgtgagggcc+ 
gcgctctccccgg+ 

20 Parameterization of the Algorithm of the Invention 

Caution: RNA/RNA and DNA/RNA duplexes contain motifs for 
wliich no literature data are available. In these cases, DNA/DNA parameters are 
assumed. Therefore, predictions might be maccurate. Users are encouraged to use 
this program with caution and discernment. 

25 No data are available for the following motifs: 

RNA/RNA single mismatches 
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RNA/DNA single mismatches 
dangling ends 
terminal mismatches 
Single mismatches 

5 Double mismatch parameters are estimated for all types of duplexes 

(DNA/DNA, RNA/RNA, DNA/RNA). 





DNA THERMODYNAMIC PARAMETERS 




Watson-Crick nearest-neishbors 






10 


SantaLucia, AUawi, and Seneviratne 

^iyyo^ rSlOCnCTnlSlry dZ>, oDDj, 

AUawi and SantaLucia (1997) 
iSiocnetnistry jo, iudoI 


12 parameters 


108 sequences 




Single mismatch nearest-neishbors 






15 


AUawi and SantaLucia (1997) 
Biochemistry 36, 100581; 
AUawi and SantaLucia (1997) 
Nucleic Acids Res., 26, 2694; Peyret 
et al., (1999) Biochemistry 38, 3468 


44 parameters 180 sequences 

Allawi and SantaLucia (1998) Biochemistry 37, 2170 
AUawi and SantaLucia (1998) Biochemistry 37, 9435 


20 


Terminal Mismatch nearest- 
neishbors 

(Appendix) 


48 parameters 


48 sequences 




Danslins end nearest-neiehbors 






25 


S. Bommarito, Peyret, SantaLucia 
(2000) Nucleic Acids Res. 28, 1929- 
1934. 


16 parameters 


16 sequences 




Na*** dependence law 








SantaLucia (1998) Proc. Natl, Acad. 
Set. USA 95, 1460-1465 


1 parameter 


86 sequences 
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RNA THERMODYNAMIC PARAMETERS 


Watson-Crick nearest-neiehbors 

Xia et al. (1998) Biochemistry 37, 14719 


12 parameters 


Sinele Mismatches 

Kierzek, Burkard, and Turner (1999) 
Biochemistry 38, 14214-14223 


44 parameters 


Terminal mismatcli nearest-neiehbors 

Freier et al. (1986) Proc. Natl Acad, Sci. USA 83, 
9373 


48 parameters 


Dangling end nearest-neishbors 

Freier et al. (1986) Proc, Natl. Acad. Sci. USA 83, 
9373 


32 parameters 


Coaxial stacking nearest-neighbors 

Walter and Turner (1994) Biochemistry 33, 12715 


16 parameters 


Loop parameters 

Matthews et al. (1999) /. MoL Biol. 288, 911-940 





HYBRID DNA/RNA THERMODYNAMIC PARAMETERS 


Watson-Crick nearest-neighbors 

Sugimoto et al. (1995) Biochemistry 34, 11211 


17 parameters 


rU»dG and rG-dT mismatches 

Sugimoto et al. (1997) Nucleic Acids Symp. Ser. 37, 199 
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DNA LOOP THERMODYNAMIC PARAMETERS 



Hairpins 

Hilbers et al. (1985) Biochie 67, 685-695 

Blommers et aL (1989) Biochemistry 28, 7491-7498 

Antao et al. (1991) Nucleic Acids Res. 19, 5901-5905 

Antaoetal. {1991) Nucleic Acids Res, 20, 819-824 

Senior et al. (1988) Proc. Natl Acad, Sci. USA 85, 6242-6246 



10 



Bulges 

LeBlanc and Morden (1991) Biochemistry 30, 4042-4047 

Zieba et al. (1991) Biochemistry 30, 8018-8026 

Ke et al. (1995) Biochemistry 34, 4593-4600 

Turner, D.H. (1992) Curr. Opin, Struc. BioL 2, 334-337 



15 



Multibranched Loops 

Kadrmas et al. (1995) Nucleic Acids Res, 23, 2122 
Lilley and Hallam (1984) /. MoL BioL 180, 179-200 

Lu et al. (1991) /. MoL BioL 223, 781-789 
Ladbury et aL (1994) Biochemistry 33, 6828-6833 
Leontis et al. (1991) Nucleic Acids Res. 19, 759-766 



20 



25 



The parameters for multibranched loops are from a best fit analysis of 
secondary structure predictions vs. experiments as done by Jaeger et aL for 
RNA (Jaeger et aL (1989) PNAS 86, 7706-7710). The current parameters for 
multibranched loops neglect the sequence and complicated length dependence 
described by Leontis and coworkers, but approximate 4-way junctions fairly 
well. Implementation of more complicated rules will require modification of 
the MFOLD algorithm. 



While the best mode for carrying out the invention has been 
described in detail, those familiar with the art to which this invention relates will 
recognize various alternative designs and embodiments for practicing the invention 
as defined by the following clamis. 
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Table 1 : Thermodynamic Parameters for Duplex Formation in IM NaCl° 



AT G AG C T C A A 
AACTCGAGTA 

AAGAGC TCTA 
ATCTCGAGAA 

AGTAGCTACA 
ACATCAATGA 

ACGATATCGA 
AGCTAT AGCA 

£TGAGCTCA£ 
£ACTCGAGT£ 

£AGAGC TCT£ 
£TCTCGAGA£ 

CGTAGC TAC£ 
£CATCGATG£ 

£CGATATCG£ 
£GCTAT AGC£ 

GTGAGCTCAfi 
fiACTCGAGTfi 

£AGAGC TCT£ 
fiT CTCGAGAfi 

fiGTAGC TACfi 
fiCATCGATGfi 

fiCGATATCGfi 
eCCTAT AGCfi 



AH"" 
(kcal / mol) 

-56.7 ± 2.5 



-55.1 ± 1.4 



-60.3 ± i.8 



-67.8 ± 1.3 



-50.6 ± 1.2 



-51.4 ± 1.3 



-55.8 ± 1.8 



-59.8 ± 1.2 



-52.6 ± 1.3 



-52.3 ± 1.4 



-59.2 ± 1.0 



-65.7 ± 0.9 



AS"*" 
(cal / mol K) 

-155.8 ± 3.7 
-149.2 ± 3.8 



-154.6 ± 4.8 



-166.9 ± 3.0 



AG 3- 
(kcal / mol) 

-9.07 ±0.12 



-8.91 ± 0.12 



-165.8 ±5.1 -9.03 ±0.12 



-192.0 ± 3.5 -8.87 ± 0.06 



-136.9 ±2.8 -8.15 ±0.11 



-137.5 ± 3.2 -8.34 ± 0.14 



-8.04 ±0.16 



-8.12 ± 0.07 



■140.7 ± 3.3 -8.57 ± 0.13 



-142.5 ± 4.1 -8.34 ± 0.08 



-164.2 ± 2.8 -8.68 ± 0.07 



-183.6 ± 2.2 



-8.81 ± 0.07 



CO 
52.3 

56.0 

54.5 

49.4 

52.9 

56.9 

49.7 

49.7 

57.7 

52.3 

51.2 

52.2 
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Table 1: Continued^ 



IT GAGC TCAl 
lACTCGAGTI 


(kcal / moJ) 
-55.4 ± 1.1 


AS"' 
(cai / mol K.) 

-149.4 ± 3.0 


AG^- ' 
(^Kcai ; nioi } 

-8.62 ±0.10 


( ^) 
56.9 


lAGAGCTCTX 
ITCTCGAGAI 


-56.5 ± 1.3 


-154.3 ± 3.7 


-8.72 ± 0.09 


54.3 


IGTAGC TACT 
ICATCGATGX 


-63.8 ± 0.8 


-178.0 ± 2.1 


-8.75 ± 0.06 


52.1 


ICGATATCGI 
TGCTAT AGCX 


-66.8 ± 0.6 


-187.9 ± 1.6 


-8.60 ± 0.04 


50.5 


£TGAGCTCAA 
AACTCGAGT£ 


-53.6 ± 1.3 


-145.8 ± 4.0 


-8.42 ± 0.06 


53.6 . 


ATGAGCTCAC 
CACTCGAGTA 


-54.0 ± 1.3 


-144.4 ± 3.2 


-8.92 ±0.15 


58.8 


£GTAGCTACA 
ACATCGATG£ 


-56.8 ± 1.4 


-155.6 ± 3.5 


-8.53 ± 0.13 


53.7 


AGTAGCTACC 
£CATCGATGA 


-57.1 ± 1.2 


-156.0 ± 3.0 


-8.71 ± 0.16 


54.4 


£CGATATCGA 
AGCTAT AGC£ 


-61.8 ± 0.5 


-172.7 ± 1.4 


-8.30 ± 0.03 


50.6 


ACGATATCG£ 
£GCTAT AGCA 


-58.3 ± 1.7 


-158.5 ± 4.2 


-8.91 ± 0.17 


56.8 


£AGAGCTCTA 
ATCTCGAGA£ 


-54.6 ± 0.6 


-147.9 ± 1.4 


-8.66 ± 0.07 


55.4 
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Table 1 : Continued^ 



AAGAGCTCT£ 
QT CTCGAGAA 



o b 



AH' 
(kcal / mol) 

-55.5 ± 1.4 



AS 



o b 



(cal / mol K) (kcal / mol) CO 
-150.4 ± 4.4 -8.85 ± 0.08 55.8 



IT GAGC TCAC 
£ACTCGAGTT 



-52.2 ± 1.0 



-140.8 ± 2.4 -8.38 ±0.11 54.9 



£TGAGCTCA1 
XACTCGAGT£ 



-55.1 ± 0.8 



-150.5 ± 1.9 -8.42 ± 0.08 53.2 



1GTAGCTAC£ -58.0 ± 1.4 
£CATCGATGT 

£GTAGC TACT -59.4 ± 0.9 
TCATCGATGC 

TC GATATCG£ -61.7 ± 1.1 
£GCTAT AGCT 

CCGATATCGI -57.9 ± 1.1 
IGCTATAGCC 

lAGAGCTCTC -55.0 ± 1.4 
CT CTCGAGAT 

CAGAGCTCTl -51.5 ± 1.1 
IT CTCGAGAC 

GT GAGC TCAA -54.2 ± 1.4 
AACGCGAGTfi 

AT GAGC TCAfi -55.4 ± 1.2 
fiACTCGAGTA 



-159.6 ±3.7 -8.39 ± 0.10 52.9 



-164.6 ± 2.1 -8.21 ± 0.06 51.7 



-170.6 ± 3.2 -8.33 ± 0.12 53.7 



-159.5 ± 2.7 -8.11 ± 0.12 52.7 



-148.0 ±3.5 -8.80 ± 0.12 57.5 



■137.9 ± 2.6 -8.44 ± 0.15 56.4 



-146.1 ± 3.3 -8.77 ± 0.16 56.8 



-148.6 ± 3.1 -9.03 ± 0.14 59.1 
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Table 1 : Continued^ 























AG°37 " 




















(kcal / mol) 


(cal / mol K) 


(kcal / mol) 


(°C) 


a G 


T 


A 


G 


C 


T A 


C 


A 


-59.4 ± 1.6 


-163.3 ± 3.9 


-S.76 ± 0.15 


53.S 


AC 


A 


T 


C 


G 


A T 


G 


Q 










A G 


T 


A 


G 


C 


T A 


C 


a 


-63.7 ± 1.6 


-174.9 ± 4.1 


-9.46 ±0.18 


56.4 


GC 


A 


T 


C 


G 


A T 


G 


A 










QC 


G 


A 


T 


A 


T C 


G 


A 


-60.4 ± 0.8 


-166.6 ± 2.1 


-8.50 ± 0.09 


53.5 


AG 


C 


T 


A 


T 


A G 


C 


Q 










A C 


G 


A 


T 


A 


T C 


G 


Q 


-61.1 ± 1.3 


-167.5 ± 3.4 


-9.04 ± 0.14 


55.8 


fiG 


C 


T 


A 


T 


A G 


C 


A 










QA 


G 


A 


G 


C 


T C 


T 


A 


-54.0 ±1.1 


-144.9 ± 2.7 


-8.82 ± 0.14 


57.9 


AT 


C 


T 


C 


G 


A G 


A 


Q 










A A 


G 


A 


G 


C 


T C 


T 


Q 


-54.8 ± 1.7 


-148.0 ± 5.0 


-8.90 ± 0.09 


56.3 


Q T 


C 


T 


C 


G 


A G 


A 


A 










I T 


G 


A 


G 


C 


T C 


A 


Q 


-56.8 ± 0.7 


-155.5 ± 1.7 


-8.64 ± 0.06 


53.5 


Q A 


C 


T 


C 


G 


A G 


T 


X 










QT 


G 


A 


G 


C 


T C 


A 


X 


-57.4 ± 0.7 


-156.7 ± 2.0 


-8.80 ± 0.05 


54.7 


X A 


C 


T 


C 


G 


A G 


T 


Q 










I G 


T 


A 


G 


C 


T A 


C 


a 


-59.2 ± 1.3 


-161.3 ± 3.6 


-8.93 ± 0.08 


56.2 




A 


T 


C 


G 


A T 


G 


X 










QG 


T 


A 


G 


C 


T A 


C 


1 


-64.8 ± 2.6 


-182.3 ± 7.2 


-8.63 ±0.17 


50.0 


I C 


A 


T 


C 


G 


A T 


G 


Q 










X c 


G 


A 


T 


A 


T C 


G 


Q 


-63.3 ± 0.6 


-176.9 ± 1.6 


-8.42 ± 0.05 


51.2 


QC 


C 


T 


A 


T 


A G 


C 


X 










QC 


G 


A 


T 


A 


T C 


G 


I 


-63.7 ± 0.9 


-177.4 ±2.2 


-8.73 ± 0.12 


52.6 


X G 


C 


T 


A 


T 


A G 


C 


G 
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Table I: Continued'. 





AH"" 
(kcal / mol) 


AS"" 
(cal/molK) 


AG",,*" 
(kcal / mol) 


(°C) 


lAGAGCTCXe 
fiXC.TCGAGAX 


-57.8 ± 0.6 


-157.4 ± 1.7 


-S.95 ± 0.03 


55.6 


GAGAGCTCTl 
ITCTCGAGAfi 


-57.3 ± 1.2 


-155.5 ± 3.3 


-8.95 ±0.10 


56.4 


Core seqiiences 

CGATATCG'^ 
GCTAT AGC 


-51.9 ± 0.6 


-145.3 ± 1.4 


-6.89 ± 0.09 


44.1 


GTAGCTAC"^ 
CATCGATG 


-51.6 ± 0.6 


-143.7 ± 1.3 


-7.01 ± 0.08 


45.6 


AGAGCTCT 
T CTCGAGA 


-50.0 ± 0.7 


-136.5 ± 1.7 


-7.76 ± 0.06 


50.2 


TGAGCTCA 
ACTCGAGT 


-50.5 ± 0.5 


-137.7 ± 1.3 


-7.73 ± 0.04 


50.4 



The top strand of each duplex is represented in the 5* to 3* orientation and the bottom strand 
is shown in the 3' to 5* direction. Tenninal mismatch nearest neighbors are represented in bold. 
Mismatches are underiined. ^ AH°, AS°, and AG^37 are the error-weighted averages of the 1/Tp^ vs. 
In Ct plot and curve fit methods in Table SI. Errors reflect the precision of the data (see text). 
Tm calculated using 10"^ M total strand concentration. ^ Data from reference (19). 
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Table 2: Nearest-neighbor thermodynamic paramters of 
like-with-like base terminal mismatches in 1 M NaCl 



Dimer 


AH° " 








(kcal/mol) 


(e.u) 


(kcal/mol) 




Terminal j 


\»A Mismatches 




A^TA 


-3.1 ± 1.3 


-7.8 ± 2.0 


-0.67 ± 0.06 


TA/AA 


-2.5 ± O.S 


-6.3 ±2,1 


-0.58 ± 0.07 


CA/GA 


-4.3 ± 1.0 


-10.7 ± 2.6 


-1.01 ± 0.07 


GA/CA 


-8.0 ± 0.7 


-22.5 ± 1.9 


-0.99 ± 0.06 




Terminal i 


Mismatches 




AC/TC 


-0.1 ± 0.6 


0.5 ± 1.5 


-0.21 ± 0.06 


T0A£ 


-0.7 ± 0.7 


-1.3 ± 1.8 


-0.29 ± 0.07 


CQ/GQ 


-2.1 ± 0.9 


-5.1 ±2.5 


-0.52 ± 0.09 


GQ/CQ 


-3.9 ± 0.7 


-10.6 ± 1.7 


-0.62 ± 0.06 




Terminal < 


G*G Mismatches 




AQ/TQ 


-1.1 ± 0.7 


-2.1 ± 1.8 


-0.42 ± 0.07 


TQ/Afi 


-1.1 ± 0.8 


-2.7 ± 2.2 


-0.29 ± 0.05 


CQ/GQ 


-3.8 ± 0.6 


-9.5 ± 1.5 


-0.83 ± 0.05 




n '7 4- A < 




— Vl.^vJ -1- VJ . \J\J 




Terminal ' 


r*T Mismatches 




Axn 


-2.4 ± 0.6 


-6.5 ± 1.6 


-0.45 ± 0.05 


TI/AI 


-3.2 ± 0.7 


-8.9 ±2.1 


-0.48 ± 0.05 


CJ/GI 


-6.1 ± 0.5 


-16.9 ± 1.2 


-0.87 ± 0.05 


Gl/CI 


-7.4 ± 0.4 


-21.2 ± LI 


-0.86 ± 0.05 



^ Thermodynamic parameters and their corresponding errors are 
calculated from Table 1 using equations 4 and 5. 
^ Dimers are given in antiparallel orientation (e.g. AC/TA equals 
5*-AC-3' paired with 3'-TA-5^). Mismatches are imderlined. 
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Table 3: Nearest-neighbor thermodynamic pWrnri^t^r^m'- 
mixed-base terminal mismatches in 1 M NaCl 

Dimer AlF^ AS^ ^ AG^T^ 

sequence^ (kcal/mol) (e.u) (kcal/mol) 



Terminal A*C Mismatches 



A A/TP 


-1.6 ± 0.7 


-4.0 ±2.1 


-0 35 ± 0 04 


AQ/TA 


-1.8 ±0.7 


-3.8 ± 1.7 


-0 59 i 0 08 


CA/G£ 


-2.6 ± 0.8 


-5.9 ± 1.8 


-0 76 i 0 07 


CC/GA 


-2.7 ± 0.7 


-6.0 ±1.6 


-0 85 ± 0 09 




-5.0 ± 0.4 


-13.8 ± 1.0 


-0 71 ±0 05 




-3.2 ± 0.9 


-7.1 ±2.2 


-1 01 ±0 10 


TA/AP 


-2.3 ± 0.5 


-5.9 ± 1.1 


-0 45 ±0 05 


TC/AA 


-2.7 ± 0.8 


-7.0 ± 2.4 


-0 55 zfc: 0 05 




Terminal C»T Mismatches 




AC/TT 


-0.9 ± 0.5 


-1.7 ± 1.4 


-0 33 ± 0 06 


AT/TC 


-2.3 ± 0.5 


-6.3 ± 1.2 


-0 35 ± 0 05 


C£/GT 


-3.2 ± 0.8 


-8.0 ± 2.0 


-0 69 i 0 07 


d/GC 


-3.9 ± 0.6 


-10.6 ± 1.2 


—\J,\J\J -E \J,\JJ 


G£/CT 


-4.9 ± 0.6 


-13.5 ± 1.7 


-0 77 -t 0 OR 


GT/CC 


-3.0 ± 0.6 


-7.8 ± 1.5 


-0 61 ±0 OR 


T£/AX 


-2.5 ± 0.8 


-6.3 ± 2.0 


-O 57 -1- 0 07 


TT/AC 


-0.7 ± 0.6 


-1.2 ± 1.6 






Terminal G«A Mismatches 






-1.9 ± 0.7 


-4.4 ± 1.8 


-0.52 ± 0.08 


Afi/TA 


-2.5 ± 0.7 


-5.9 ± 1.7 


-0.65 ± 0.07 




-3.9 ± 0.8 


-9.6 ± 2.1 


-0.88 ± 0.09 


CQ/GA 


-6.0 ± 0.9 


-15.5 ± 2.1 


-1.23 ±0.1 




-4.3 ± 0.5 


-11.1 ± 1.3 


-0.80 ± 0.06 


GQ/CA 


-4.6 ± 0.7 


-11.4 ± 1.8 


-1.08 ± 0.09 


TA/AQ 


-2.0 ± 0.7 


-4.7 ± 1.6 


-0.53 ± 0.07 


Tfi/AA 


-2.4 ± 0.9 


-5.8 ± 2.7 


-0.57 ± 0.05 
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Table 3: Continued 


Dimer 


AH"" 


AS"" 


AG°37'' 


sequence^ 


(kcal/mol) 


(e.u) 


(kcal/mol) 




Terminal G»T Mismatches 




Afi/TT 


-3.2 ± 0.4 


-8.7 ± 1.1 


-0.45 ± 0.04 


AX/TQ 


-3.5 ± 0.4 


-9.4 ± 1.2 


-0.54 ± 0.03 


Cfi/GI 


-3.8 ± 0.7 


-9.0 ± 1.9 


-0.96 ± 0.06 


CI/GO 


-6.6 ± 1.3 


-18.7 ±3.6 


-0.81 ± 0.09 


GQ/CX 


-5.7 ± 0.4 


-15.9 ± 1.0 


-0.76 ± 0.05 


GT/CQ 


-5.9 ± 0.5 


-16.1 ± 1.3 


-0.92 ± 0.07 


TQ/AX 


-3.9 ± 0.5 


-10.5 ± 1.2 


-0.59 ± 0.03 


TX/AQ 


-3.6 ± 0.7 


-9.8 ± 1.9 


-0.59 ± 0.06 



^ Thennodynamic parameters and their corresponding errors 
are calculated from Table 1 using equations 4 and 5. 

^ Dimers are given in antiparallel orientation (e.g. AQ/TA equals 
5'-AC-3' paired with 3 -TA-5*). Mismatches are underlined. 
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Table S 1 : Thermodynamic Parameters for Duplex Formation in IM NaCl^ 

























AS' 




AG" 


1 
























(kcal / mol) 


(cal / mol K) 


(kcal / mol) 


(°C) 


A T 






VJ 


c 


T 


c 


A 


A 


c 


-56.0 ± 


3.2 


-151.4 ± 


10.1 


-9.07 ± 


0.13 


57.0 


A A 


C 


T 


c 


G 


A 


G 


T 


A 


d 


-57.6 ± 


4.0 


-156.5 ± 


4.0 


-9.11 ± 


0.46 


56.7 


A A 


G 


A 


G 


c 


T 


c 


T 


A 


c 


-54.6 ± 


1.8 


-146.9 ± 


5.3 


-8.91 ± 


0.12 


57.7 


AT 


C 


T 


C 


G 


A 


G 


A 


A 


d 


-56.0 ± 


2.3 


-151.6 ± 


5.4 


-8.97 ± 


0.64 


56.4 


A G 


T 


A 


G 


c 


T 


A 


c 


A 


c 


-59.9 ± 


2.2 


-164.0 ± 


6.6 


-9.03 ± 


0.12 


55.4 


A C 


A 


T 


C 


A 


A 


T 


G 


A 


d 


-61.3 ± 


3.3 


-168.3 ± 


7.9 


-9.09 ± 


0.83 


55.3 


A C 


G 


A 


T 


A 


T 


c 


G 


A 


c 


-60.5 ± 


2.4 


-166.5 ± 


7.7 


-8.86 ± 


0.06 


54.3 


A G 


C 


T 


A 


T 


A 


G 


C 




d 


-70.8 ± 


1.5 


-198.4 ± 


3.9 


-9.28 ± 


0.34 


53.6 


C T 


G 


A 


G 


c 


T 


c 


A 


r 


c 


-52.3 ± 


4.0 


-142.4 ± 


12.7 


-8.16 ± 


0.11 


52.4 


£ A 


C 


T 


C 


G 


A 


G 


T 


£ 


d 


-50.4 ± 


1.2 


-136.6 ± 


2.8 


-8.06 ± 


0.36 


52.4 


C A 


G 


A 


G 


C 


T 


c 


T 




c 


-55.6 ± 


2.5 


-152.1 ± 


7.7 


-8.36 ± 


0.14 


52.8 


£ T 


C 


T 


C 


G 


A 


G 


A 


£ 


d 


-49.9 ± 


1.5 


-134.5 ± 


3.5 


-8.14 ± 


0.45 


53.1 


£ G 


T 


A 


G 


C 


T 


A 


C 


£ 


c 


-55.1 ± 


2.4 


-151.8 ± 


7.2 


-8,04 ± 


0.16 


50.9 


£ C 


A 


T 


C 


G 


A 


T 


G 


£ 


d 


-56.8 ± 


2.7 


-156.9 ± 


6.6 


-8.08 ± 


0.70 


50.7 


£ C 


G 


A 


T 


A 


T 


c 


G 


c 


c 


-57.3 ± 


2.6 


-158.4 ^ 


8.1 


-8.12 ± 


0.07 


50.8 


£ G 


C 


T 


A 


T 


A 


G 


C 


£ 


d 


-60.5 ± 


1.3 


-168.3 ± 


3.3 


-8.21 ± 


0.31 


50.9 


G T 


G 


A 


G 


c 


T 


c 


A 


ri 


c 


-55.2 ± 


2.1 


-150.4 ± 


6.4 


-8.58 ± 


0.14 




QA 


C 


T 


C 


G 


A 


G 


T 


£ 


d 


-50.9 ± 


1.7 


-137.1 ± 


3.9 


-8.38 ± 


0.48 


54.4 


£ A 


G 


A 


G 


C 


T 


C 


T 


£ 


c 


-56.8 ± 


2.0 


-155.1 ± 


5.6 


-8.72 ± 


0.26 


54.6 


Q T 


C 


T 


C 


G 


A 


G 


A 


£ 


d 


-48.0 ± 


1.9 


-128.1 ± 


5.9 


-8.30 ± 


0.09 


54.9 



55 



wo 01/94611 



PCT/USOl/18424 



Table SI: Continued. 



QG 7 A G C T AC G. 
GCATCGATGfi 

GCGATATCGG 
QG C T A T AG C G 



AH° 
(kcal / mol) 

-56.9 ± 1.5 



-61.1 



1.4 



-62.4 ± 2.4 
-66.2 ± 0.9 



(cal / mol K) 

-155.4 ± 4.8 
-168.5 ± 3.4 



AG, J 
(kcal / mol) 

-8.67 ± 0.07 
-8.83 ± 0.35 



-172.9 ± 7.4 -8.79 ± 0.08 
-184.7 ± 2.3 -8.93 ± 0.21 



(°C) 

54.2 
53.9 

53.3 
53.0 



ITGAGCTCAI 
lACTCGAGTl 

TAGAGCTCTT 
ITCTCGAGAX 

IGTAGCTACl 
TCATCGATGX 

ICGATATCGX 
JGCTAT AGCl 

£TGAGCTCAA 
AACTCGAGT£ 

ATGAGCTCAC 
CACTCGAGTA 

£GTAGC TACA 
ACATCGATGC 

AGTAGCTACC 
CCATCGATGA 

CCGATATCGA 
AGCTAT AGC£ 



" -56.8 ± 1.5 

-53.7 ± 1.7 

' -58.7 ± 1.6 

^ -52.9 ± 2.1 

' -63.2 ± 1.1 

-64.4 ± 1.1 

= -64.6 ± 1.3 

-67.4 ± 0.7 

' -56.0 ± 2.1 

-52.0 ± 1.7 

" -57.0 ± 2.4 

^ -52.7 ± 1.6 

" -57.8 ± 2.9 

" -56.5 ± 1.6 



-59.4 ± 4.0 
-56.8 ± 1.3 



-61.6 ± 1.1 
-61.9 ± 0.6 



-155.1 ± 4.7 

-145.6 ± 3.9 

-160.3 ± 4.6 

-142.5 ± 6.4 

-175.5 ± 3.4 

-179.4 ± 2.6 

-180.4 ± 4.0 

-189.4 ± 1.8 

-153.1 ± 6.2 

-140.6 ± 5.3 

-154.9 ± 7.1 

-141.7 ± 3.6 

-158.8 ± 8.9 

-155.0 ± 3.8 

-163.2 ± 12.2 

-155.5 ± 3.1 

-171.8 ± 3.5 

-172.9 ± 1.5 



-8.63 ± 0.10 54.2 

-8.52 ± 0.46 54.4 

-8.96 ± 0.19 55.4 

-8.66 ± 0.10 55.6 

-8.75 ± 0.06 52.9 

-8.80 ± 0.25 52.8 

-8.59 ± 0.05 51.7 

-8.69 ± 0.16 51.5 

-8.58 ± 0.15 53.9 

-8.39 ± 0.07 54.1 

-8.94 ± 0.16 55.8 

-8.73 ± 0.46 56.1 

-8.54 ± 0.13 53.2 

-8.46 ± 0.42 53.0 

-8.74 ± 0.18 53.9 

-8.60 ± 0.35 53.8 

-8.30 ± 0.03 50.8 

-8.30 ± 0.14 50.7 
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Table SI: Continued. 















AS' 


0 


AG^ 


> 

37 


■-p b 
'■M 










(kcal / mol) 


(cal / mol K) 


(kcal / mol) 


(°C) 


A C 


G A T A 


T C G £ 


C 


-61.8 


± 


3.2 


-170.3 


± 


9.6 


-8.94 ± 


0.19 


54.3 


C G 


C T A T 


A G C A 


d 


-57.0 


db 


1.9 


-155.8 


± 


4.6 


-8.71 ± 


0.52 


54.4 


£ A 


G A G C 


T C T A 


c 


-55.7 


± 


1.1 


-151.6 


± 


3.5 


-8.67 ± 


0.07 


54.6 


A T 


C T C G 


A G A £ 


d 


-54.2 


i 


0.7 


-14/.Z 


_ r 


l.o 




U. lo 




A A 


G A G C 


T C T £ 


c 


-58.6 


± 


4.1 


-159.8 


± 


12.4 


-9.07 ± 


0.26 


56.1 


£ T 


C T C G 


A G A A 


d 


-55.1 


± 


1.5 


-149.1 


± 


4.7 


-8.83 ± 


0.08 


55.8 


I T 


G A G C 


T C A £ 




-54.1 


± 


1.9 


-147.4 




5.8 


-8.40 ± 


0.12 


53.4 


£ A 


C T C G 


A G T X 


d 


-51.5 




1.1 


-139.4 




2.6 


-8.27 ± 


0.32 


53.5 


C T 


G A G C 


T C A T 


c 


-57.3 


± 


4.5 


-157.5 




14.3 


-8.43 ± 


0.09 


52.7 


1 A 


C T C G 


A G T £ 


d 


-55.0 




0.8 


-150.3 


± 


1.9 


-8.37 ± 


0.22 


53.0 


I G 


T A G C 


T A C £ 


c 


CO n 
05./ 


ziz 


2.1 


-162.3 


ifc 


6.5 


-8.39 ± 


0.10 


52.0 


£ C 


A T C G 


A T G I 


d 


-57.4 


± 


1.9 


-158.2 


± 


4.6 


-8.32 ± 


0.48 


52.0 


£ G 


T A G C 


TACT 


c 


-59.3 


db 


1.6 


-160.7 




5.1 


-8.21 ± 


0.07 


58.1 


1 C 


A T C G 


A T G £ 


d 


-59.5 


± 


1.2 


-165.4 


± 


2.3 


-8.19 ± 


0.28 


50.7 


1 C 


G A T A 


T C G £ 


c 


^fkO ft 

-OZ..O 


4- 




1 7^ 7 


,1. 








50.7 


£G 


C T A T 


A G C X 


d 


-57.7 


± 


2.4 


-159.6 


± 


5.7 


-8.16 =t 


0.59 


51.0 


£ C 


G A T A 


T C G I 


c 


-61.1 


± 


1.8 


-170.6 


± 


5.5 


-8.13 ± 


0.13 


50.0 


1 G 


C T A T 


A G C £ 


d 


-56.3 


zb 


1.3 


-155,8 




3.1 


-7.98 ± 


0.34 


50.2 


I A 


G A G C 


T C T £ 


c 


-56.4 


± 


2.0 


-153.5 


± 


6.3 


-8.81 ± 


0.12 


55.2 


£ T 


C T C G 


A G A I 


d 


-53.8 


± 


1.8 


-145.5 




4.3 


-8.69 db 


0.52 


55.4 


£ A 


G A G C 


T C T I 


c 


-54.8 


± 


2.1 


-149.4 


± 


6.4 


-8.47 ± 


0.17 


53.7 


1 T 


C T C G 


A G A £ 


d 


-50.3 


± 


1.2 


-135.6 


± 


2.8 


-8.28 ± 


0.36 


53.9 
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Table SI: Continued. 



fiTGAGCTCAA 
AACGCGAGTS 

ATGAGCTCAfi 
SACTCGAGTA 

fiG TAGC 7 AC A 
ACATCGATGG 

AGTAGC TACfi 
fiCATCGATGA 

fiCGATATCGA 
AGCTAT AGCe 

ACGATATCGfi 
fiGCTAT AGCA 

fiAGAGCTCTA 
AT CTCGAGAS 

AAGAGCTCTS 
fiTCTCGAGAA 

IT GAGCTCAfi 
fiACTCGAGTl 

GT GAGC TCAI 
XACTCGAGTfi 

IGTAGCTACfi 
fiCATCGATGX 

SGTAGCTACI 
XCATCGATGS 



(kcal / mol) 

" -56.6 ± 2.8 

-53.5 ± 1.6 

^ -57.7 ± 1.9 

^ -53.8 ± 1.6 

' -59.2 ± 4.2 

-59.4 ± 1.7 

"= -63.3 ± 3.4 

^ -63.8 ± 1.8 

' -62.1 ± 1.3 

^ -59.1 ± 1.1 

" -63.1 ± 2.6 

-60.4 ± 1.6 

" -56.4 ± 1.9 

-52.8 ± 1.3 

" -57.1 ± 2.6 

^ -53.1 ± 2.2 

" -55.2 ± 1.5 

-57.2 ± 0.8 

' -57.3 ± 0.7 

-57.6 ± 1.6 



-59.9 ± 1.8 
-58.2 ± 2.0 



-62.8 ± 3.6 
-67.1 ± 3.8 



AS° 
(cal / mol K) 

-154.3 ± 8.4 
-144.6 ± 3.6 

-157.0 ± 5.6 
-145.1 ± 3.6 

-162.7 ± 12.9 
-163.3 ± 4.1 

-173.4 ± 10.2 
-175.2 ± 4.5 

-172.8 ± 3.8 
-163.6 ± 2.6 

-174.1 ± 7.9 
-166.0 ± 3.7 

-153.5 ± 5.6 
-142.4 ± 3.1 

-154.7 ± 7.6 
-142.8 ± 6.7 

-150.0 ± 4.5 
-156.4 ± 1.8 

-156.5 ± 2.3 
-157.4 ± 3.7 

-164.5 ± 5.5 
-158.9 ± 4.8 

-174.8 ±11.0 
-187.8 ± 9.4 



AG°37 

(kcal / mol) (°C) 

-8.79 ± 0.17 55.1 

-8.64 ± 0.44 55.2 

-9.05 ± 0.15 56.3 

-8.85 ± 0.45 56.5 

-8.77 ± 0.16 54.1 

-8.72 ± 0.44 53.8 

-9.46 ± 0.20 56.8 

-9.47 ± 0.45 56.6 

-8.51 ± 0.10 51.8 

-8.39 ± 0.27 51.9 

-9.06 ± 0.15 54.6 

-8.92 ± 0.40 54.6 

-8.84 ± 0.15 55.4 

-8.67 ± 0.38 55.6 

-9.08 ± 0.20 56.7 

-8.86 ± 0.10 56.8 

-8.63 ± 0.07 54.5 

-8.71 ± 0.21 54.3 

-8.80 ± 0.05 54.9 

-8.81 ± 0.42 54.8 

-8.93 ± 0.08 54.8 

-8.86 ± 0.53 55.0 

-8.62 ± 0.17 52.3 

-8.79 ± 0.85 52.1 
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Table SI: Continued. 



ICGATATCGfi 
fiCCTATAGCl 

fiCGATATCGl 
TGCTAT AGCfi 



I A G A G C 



T n T cz 



QT CTCGAGAX 

CAGAGCTCTI 
IT CTCGAGAS 



AH° 
(kcal / mol) 

" -58.9 ± 2.8 
^ -63.5 ± 0.6 

" -61.6 ± 3.7 
^ -63.9 ± 0.9 

' -5^.4 ± 1.5 
-58.0 ± 0.6 

" -57.7 ± 1.6 
-56.6 ± 1.9 



AS° 
(cal / mol K) 

-163.0 ± 8.9 
-177.3 ± 1.6 

-170.4 ± 11.5 
-177.6 ± 2.3 

-153.3 ± 4.5 
-158.1 ± 1.9 

-157.2 ± 4.7 
-153.9 ± 4.6 



AG°37 

(kcal / mol) 

-8.40 ± 0.06 
-8.53 ±0.15 

-8.71 ± 0.14 
-8.77 ± 0.22 

-8.89 ± 0.06 
-8.96 ± 0.03 

-8.95 ±0.10 
-8.90 ± 0.52 



CO 

52.0 
51.6 

53.1 
52.8 

55.7 



55.7 
55.7 



Core sequences 

CGATATCGe = .55.7 ± 3.9 -157.1 ± 12.1 -6.93 ± 0.12 44.1 

GCTATAGC -51.8 ± 0.6 -145.1 ± 1.4 -6.82 ± 0.15 44.0 

GTAGCTAC" " -55.1 ± 2.3 -155.0 ± 7.0 -7.04 ± 0.10 44.9 

CATCGATG -51.4 ± 0.6 -143.3 ± 1.3 -6.95 ± 0.14 44.9 

AGAGCTCT " -49.5 ± 1.8 -134.5 ± 5.7 -7.76 ± 0.07 50.6 

TCTCGAGA ^ -50.1 ± 0.8 -136.7 ± 1.8 -7.76 ± 0.22 50.5 

TGAGCTCA " -50.7 ± 0.7 -138.4 ± 2.2 -7.73 ± 0.04 50.1 

ACTCGAGT -50.3 ± 0.7 -137.3 ± 1.6 -7,72 ± 0.18 50.1 



^ The top strand of each duplex is represented in the 5' to 3' orientation and the bottom strand 

is shown in the 3' to 5' direction. Terminal mismatch nearest neighbors are represented in bold. 
Mismatches are underlined. T^, calculated using 10"^ M total strand concentration. 

Thermodynamic parameters from averaging the fits of melting curves. Reported errors are 
standard deviations in the precision of the data. Thermodynamic parameters from T^" vs. \n(Cj) 

plots. Reported errors are standard deviations in the precision propagated from the slope and 
intercept of the l/T^ vs. In Cj plot. * Data from reference (19). 
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Table 1 : Thermodynamic Parameters for Hairpin Oligomer Association and Oligomer Duplex 
Formation. 

(kcal / mol) (cal / mol K) (kcaJ mol ) i ^'c » 

Slvqtems ^^-ith Elementary Interfaces^ 



C 



C 



AAGCCTTGA-ACAACG ' -50.7 ± 4.1 -141.7 ± 11.3 -6.74 = 0.27 4Q.3 
CGCGGAACT/TGTTGC^'* 

AAGCCTTGT-TCAACG " -52.8 ± 4.2 -146.3 ± 11.7 -7.43 a 0.30 53.1 
CGCGGAACA/AGTTGC ^"^ 



GCAACT-TGTTCCGAA^ = -63.5 ± 5.1 -179.2 ± 14.3 -7.99 ± 0.32 53.2 
C G T TG A / A C AAGGCC C-^ 



C 



C 



AAGCCTTGA-TCAACG "-53.6* 4.3 -149.3 * 11.9 -7.34 * 0.29 52.3 



(iiO 



c 



c 



CGCGGAACT/AGTTGC 

AAGCCTTGT-ACAACG '=-45.1 * 3.6 -124.8 * 10.0 -6.42 * 0.26 48.4 
CGCGGAACA/TGTTGC 



^A A GC C T T G C - AC AAC G * -46.1 * 3.7 -128.5 * 10.3 -6.26 * 0.25 46.9 
^CGCGGAACG/TGTTGC 

^AAGCCTTGT-GCAACG ' -52.2 * 4.2 -144.4 * 11.5 -7.39 * 0.30 53.1 
^CGCGGAACA/CGTTGC 

AAGCCTTGG-TCAACG * -53.6 * 4.3 -148.2 * 11.9 -7.67 * 0.31 54.4 
CGCGGAACC/AGTTOC 

A AGC C TT G A - CCAACG * -51.3 * 4.1 -140.2 * 11.2 -7.81 * 0.31 56.2 



CGCGGAACT/GGTTGC 
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Table 1 : Continued. 

AH** AS° AG°3- Ts, 

(kcal/mol) (cal/molK) (kcal mol) (^"Ci 

^ A A GC C T T G C - T C AAC G '-46.1*3.7 -126.3 ± 10.1 -6.90 = 0.28 SI." 
^CGCGGAA CG/AGTTGC 

A AGC C T T G A - GC AAC G ' -48..'2 ± 3.9 -131.7 ± 10.5 -7.32 ± 0.29 53.9 



c 



CGCGGAACT/CGTTGC 



GCAACA-GGTTCCGAA-v^ -51.9 ± 4.2 -147.0 ± 11.8 -6.34 ± 0.25 46.3 
CGTTGT/CCAAGGCCC>' 

^A A GC C T T G G - AC AAC G ' -50.4 ± 4.0 -139.5 ± 11.2 -7.10 ± 0.28 51.7 
^CGCGGAACC/TGTTGC 

^AAGCCTTGT-CCAACG * -54.2 ± 4.3 -147.9 ± 11.8 -8.29 ± 0.33 58.3 
^C GCGGAA C A / GGTTGC 

-AAGCCTTGC-GCAACG " -47.6 ± 3.8 -130.9 ± 10.5 -6.96 ± 0.28 51.6 
^CGCGGAACG/CGTTGC 

^AAGCCTTGG-CCAACG ' -52.1 ± 4.2 -140.1 ± 11.2 -8.67 ± 0.35 61.9 
GCGGA A C C / GGTTGC 

^AAGCCTTOC-CCAACG ' -53.3 ± 4.3 -146.5 * 11.7 -7.92 ± 0.32 56.1 
^CGCGGAACG/GGTTGC 

^AAGCCTTGG-GCAACG * -49.3 ± 3.9 -135.2 ± 10.8 -7.35 ± 0.29 53.8 
^ C G C G G A ACC/CGTTGC 
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Table 1: Continued. 



AH" 
(kcal / mol) 

jsysTems with Dangling Ends at TnTffrfHffff" 

^AAGCCTTGC-GCAACG ' -44.4 ± 3.6 
^CGCGGAA CG/CGTTGC 

A 

^AAGCCTTGC-GCAACG * -48.0 ± 3.8 
^CGCGGAACG/CGTTGC 



AS° AG°3"» Tsi 

(cal / moi K) (kcal /mol) (X ) 



•121.4 ± 9.7 -6.76 ± 0.27 51.1 



•131.8 ± 10.5 -7.16 ± 0.29 52.9 



c 



AAGCCTTGC-GCAACG ' -46.3 ± 3.7 -127.4 ± 10.2 -6.83 ± 0.27 51.0 



CGCGGAACG/CGTTGC 
T 



c 



AAGCCTTGC-GCAACG " -49.0 ± 3.9 -136.6 ± 10.9 -6.59 ± 0.26 48.7 



CGCGGAACG/CGTTGC 

T 



c 



AAGCCTTGC-GCAACG * -37.6 ± 3.0 -102.0 ± 8.2 -5.91 ± 0.24 46.4 



CGCGGAACG/CGTTGC 
A A 



^AAGCCTTGC-GCAACG * -36.2 ± 2.9 -97.3 ± 7.8 -6.03 ± 0.24 47.8 



CGCGGAACG/CGTTGC 
T T 



c 



AAGCCTTGC-GCAACG ' -44.0 ± 3.5 -123.5 * 9.9 -5.67 * 0.23 43.0 



CGCGGAACG/CGTTGC 
A T 



^AAGCCTTGC-GCAACG -43.6± 3.5 -119.5 ± 9.6 -6.53 ± 0.26 49.6 



CGCGGAACG/CGTTGC 
T A 



^A AGC C T T GG - TCAACG « -47.2 ± 3.8 -129.8 * 10.4 -6.90 * 0.28 51.3 
^CGCGGAACC/AGTTGC 

A 
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Table 1: Continued. 



(kcal / mol) (cal / mol K) (kcal / moi) (°C) 
Systems with Extra Central NuclentiHe at the Interface" 

AAGCCTTGCAGCAACG ' -44.4 ± 3.6 -122.2 ± 9.8 -6.50 ± 0.26 49.2 



c 



CGCGGAACG/CGTTGC 



^AAGCCTTGTACCAACG ' -45.0 ± 3.6 -124.6 ± 10.0 -6.32 ± 0.25 47.7 
^CGCGGAACA/GGTTGC 



Oligomers 



TCAACG ' -38.5 ± 3.1 -108.3 ± 8.7 -4.94 ± 0.20 38.0 
AGTTGC 



ACAACG * -36.1 ± 2.9 -99.3 ± 7.9 -5.32 ± 0.21 41.4 
TGTTG C 



GCAACG ' -42.5 ± 3.4 -117.1 ± 9.4 -6.21 ± 0.25 47.5 
CGTTG C 



CC AAC G ' -38.6 ± 3.1 -106.2 ± 8.5 -5.69 ± 0.23 44.2 
GGTTGC 



TGTTG C ' -37.1 ± 3.0 -101.1 ± 8.1 -5.79 ± 0.23 45.3 
ACAACG 



AGTTGC * -37.0 db 3.0 -lOl.O ± 8.1 -5.70 ± 0J23 44.2 
TCAAC O 
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^ Tm calculated using 1 0"^ total strand concentration. 

^The top strand of each system is conventionally represented in the 5* to 3' orientation. Nucleotides 
involved in coaxial stacking interfaces are represented in bold. 

Parameters obtained by averaging the results of melt fit and TM"* vs. ln(CT/4) plot methods. 
Errors are estimated to be 8% for AH° and AS^ and 4% for AG'^j?. 
(I). (II). (Ill) p^eiting curves for these systems are shown in Figure 2. 
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Jl 11 -^U .IJ U U 41 .."-tM 



Table 2: Thermodynamic Parameters for Coaxial Stacking'. 



AH*'(coaxial stacking) 
(kcal / mol) 



Fl>»Tnpntarv Interfaces 

G A - AC 
CT / T G 

GT - T C 
C A / AG 

CT - TG 
GA / AC 

GA - TC 
CT / AG 

GT - AC 
C A / TG 

GC - AC 
CG / TG 

GT - GC 

C A / CG 

GG - TC 
CC / AG 

GA - CC 
CT / GG 

GC - TC 
CG / AG 

GA - GC 
CT / CG 

C A - GG 
GT / CC 

GG - AC 
CC / TG 



-14.6 ± 5.0 

-14.3 ± 5.2 

-26.6 * 5.9 

-15.1 ± 5.3 

-9.0 ± 4.6 

-10.0 ± 4.7 

-9.6 ± 5.4 

-15.1 ± 5.3 

-12.7 ± 5.1 

-7.6 ± 4.8 

-5.6 ± 5.1 

-14.8 ±. 5.1 

-14.2 ± 5.0 



AS°(coaxial stacking) AG''37(coaxial stacking) 
(cal / mol K) (kcal / mol) 



-42.4 ± 13.8 

-38.0 ± 14.6 

-78.2 ± 16.5 

-41.0 ± 14.8 

-25.5 ± 12.8 

-29.2 ± 13.0 

-27J ± 14.9 

-39.9 ± 14.7 

-34.0 ± 14.1 

-18.0 ± 13.3 

-14.6 ± 14.1 

-45.9 ± 14.3 

-40.2 ± 13.7 



-1.42 ± 0.34 

-2.49 ± 0.36 

-2.29 ± 0.39 

-2.40 ± 0.35 

-1.10 ± 0.33 

-0.94 ± 0.33 

-1.18 ±. 0.39 

-2.73 ± 0.36 

-2.12 ± 0.39 

-1.97 ± 0.34 

-1.11 ± 0.38 

-0.56 ± 0.34 

-1.78 ± 0.35 
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Table 2: Continued. 



AH^Ccoaxial stacking) AS°(coaxial stacking) AG°37(coaxial stacking) 

(kcal / mol) (cal / mol K) (kcal / mol) 

GT-CC -15.6 ± 5.3 -41.8 ± 14.6 -2.61 ± 0.40 
CA / GG 

G C - G C -5.0 ± 5.1 -13.8 ± 14.1 -0.75 ± 0.37 
CG / CG 

GG- CC -13.5 ± 5.2 -33.9 ± 14.1 -2.98 ± 0.41 
CC / GG 

G C - C C -14.7 ± 5.3 -40.3 ± 14.5 -2.23 ± 0.39 
CG / GG 

GG-GC -6.8 ± 5.2 -18.1 ± 14.3 -1.14 '± 0.38 
CC / CG 

Interfaces with Dangling Ends'* 

GC-GC -1.9 ± 4.9 -4.3 ± 13.5 -0.55 ± 0.37 
CG / CG 
A 

GC-GC -5.5 ± 5.1 -14.7 ± 14.1 -0.95 ± 0.38 
CG / CG 
A 

GC-GC -3.8 ± 5.0 -10.3 ± 13.8 -0.62 ± 0.37 
CG / CG 
T 

GC-GC -6.4 ± 5.2 -19.5 ± 14.4 -0.38 ± 0.36 
CG / CG 
T 

GC-GC 5.0 ± 4.5 15.1 ± 12.4 0.30 ± 0.34 
CG / CG 
A A 
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Table 2: Continued. 



AH°(coaxiaJ stacking) AS°(coaxial stacking) AG°37(coaxial stacking) 

(kcal/mol) (cal/molK) (kcal/mol) 

GC-GC 6.3 ± 4.5 19.8 ± 12.2 0.18 ± 0.35 

CG / CO 
T T 

GC-GC -1.4 .± 4.9 -6.4 ± 13.6 0.55 ± 0.34 

CG / CG 
A T 

GC-GC -1.1 ± 4.9 -2.4 ± 13.4 -0.32 ± 0.36 

CG/CG 

T A 

A 

GG-TC -8.6 ± 4.9 -21.5 ± 13.5 -1.96 ± 0.34 

CC / AG 
A 

Interface with Extra Central NuelentiHe* 

GCAGC -1.9 ± 4.9 -5.1 ± 13.5 -0.29 ± 0.36 

CG/CG 

GTACC -6.4 ± 4.7 -18.5 ± 13.1 -0.64 ± 0.34 

CA / GG 



* These parameters and their coiresponding errors are deduced firom Table 1 as described in the text. 
** The top strand of each duplex is conventionally represented in the 5' to 3' orientation. Nucleotides 
involved in coaxial stacking interfaces are represented in bold. 
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Table SI: Extinction coefficients of hairpins at 25 °C 



^AAGCCTTGGTCAACG 
GCGGAACC 



(L mor* cm**) (L mof cm'') 
188847 192310 



calculated 



C 



AAGCCTTGATCAACG 
C GCGGA ACT 



188718 



195950 



c 



AAGCCTTGCTCAA CG 
CGCGGAACG 



186139 



191610 



C 



AAGCCTTGTTCAACG 
C GCGG A ACA 



191071 



195810 



C 



AAGCCTTGAACAACG 
C GCGGA ACT 



194330 



200950 



C 



AAGCCTTGTACAACG 
C GCGGA ACA 



193953 



200S10 



C 



AAGCCTTGGACAACG 
C GCGGAACC 



188889 



197310 



C 



AAGCCTTGGACAACG 
CGCGGAACG 



194623 



196610 



C 



A AGC C T TGAGC AACG 
C GCGGAACT 



1929S3 



197350 



c 



AAGCCTTGTGCAACQ 
C GCGGAACA 



195968 



197630 



C 



AAGCCTTGGGCAACG 
CGCGGAACG 



195177 



193710 
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Table SI: Continued 



AAGCCTTGCGCAACG 
•^CGCGGAACG 



190944 



193010 



C 



AAGCCTTGACCAACG 
C GCGGA ACT 



192663 



194350 



C 



AAGCCTTGTCCAACG 
C GCGGA A C A 



193532 



194210 



^AAGCCTTGGCCAACG 
^C GCGG AACC 



195094 



192390 



C 



AAGCCTTGCCCAACG 
C GCGGA ACG 



192864 



190010 



GCAACA-G TTCC A A 
CC AAGGCC C 



200806 



196010 



GCAACT-T TTCC A A 
AC AAGGCC C 



206650 



192810 



c 



AAGCCTTGCAGCAACG 
C GCGGA ACG 



202688 



206510 



^AAGCCTTGTACCAACG 
^C GCGGA ACA 



204450 



208010 



^AAGCCTTGCGCAACG 
^C GCGGAACG 



191282 



193010 



C 



AAGCCTTGCGCAACG 
C GCGGAACG 
T 



198343 



193010 
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^ Calculated wth Equation 3. 
^ Calculated with Equation 4. 
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Table S2: Thermodynamic Parameters for Hairpin Oligomer Association and Oligomer Dupiex 
Formation. 

AH° AS** AG°3: Ts,' 

(kcal/mol) (cal/mol-K) (kcal/mol) ("^C) 

Flemgntar^- inrerfages' 



C 



A A GC C T T G A - AC AA C G (i) c ^8.1 ± 2.5 -133.2 ± 8.1 -6.8^0.1 50.1 



CGCGGAACT/TGTTGC "-53.3*2.7 -150.1 ± 8.7 -6.7* 0.0 48.4 



^AAGCCTTGT - TCAACG («) c -55.7 ± 4.3 -155.8 ± 14.0 -7.4 ± 0.2 52.1 

^CGCGGAACA/AGTTGC "-49.9* 1.7 -136.9 ± 5.6 -7.4 ± 0.0 54.2 

^A AGC C T T G A - TC AA C G(«0 e -51.9 ± 4.0 -143.6 ± 12.7 -7.3 ± 0.1 52.8 

^CGCGGAACT/AGTTGC "-55.4* 2.1 -155.0 * 4.9 -7.3 * 0.5 51.7 

GCAACT-TGTTCCGAA-s « -61.6 * 2.6 -173.1 * 7.9 -8.0 * 0.1 53.7 

CGTTGA/ACAAGGCCC^ "-65.5 * 3.1 -185.2 * 9.8 -8.0 * 0.0 52.8 



C 



AAGCCTTGT-ACAACG e -48.0 * 2.3 -134.4 * 7.2 -6.4 * 0.1 47.3 



CGCGGAACA/TGTTGC "-42.2* 1.9 -115.2* 6.3 -6.5 *0.l 49.5 



-AAGCCTTGC-ACAACG c ^5.9 ± 1.1 -127.9 * 3.3 -6.3 * 0.1 47.1 

^CGCGGAACG/TGTTGC " -46.3 * 4.0 -129.2 * 13.1 -6.3 * 0.1 46.8 

-AAGCCTTGT-GCAACG «= -54.7 * 2.7 -152.4 * 9.3 -7.4 * 0.2 52.2 

^CGCGGAACA/CGTTGC ** -49.7 * 3.4 -136.3 * 11.0 -7.4 * 0.1 53.9 



AAGCCTTGG-TCAACG » -53.5 * 4.3 -147.7 * 13.6 -7.7 * 0.2 54.4 
CGCGGAACC/AGTTGC " -53.8 * 3.4 -148.7 * 11.1 -7.7 * 0.1 54.3 



^A AGC C T T G A - CC A A C G c .52.9 * 2.2 -145.4 * 7.3 -7.8 * 0.1 55.6 
^CGCGGAACT/GGTTGC " -49.6 * 2.4 -134.9 * 5.5 -7.8 * 0.7 56.8 



A AG CCTTGC-TCAACG « -48.6 * 2.8 -134.6 * 9.3 -6.9 * 0.2 50.6 
CGCGGAACG/AGTTGC "-43.6 * 2.2 -118.1 * 7.3 -6.9 * 0.1 52.8 



-A AGC C T T GA - GC A A C G c -49.7 ± 3.5 -136.6 * 11.4 -7.3 * 0.1 53.4 
^CGCGGAACT/CGTTGC " -46.6 * 1.8 -128.8 * 5.7 -7.3 * 0.0 49.9 
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Table S2: Continued. 



(kcal / mol) 

GCAACA-GGTTCCGAA-. ' -53.1 x 3.8 
CGTTGT/CCAAGGCCC-^ " -50.8 ± 3.3 



AS" 
(cal / mol K) 

■ 150.8 3: 13.0 
-143.3 ± 10.9 



(kcal ' mol) ^ 



-6.3 = 0.2 
-6.4 = 0.1 



45.8 
46.9 



c 



AAGCCTTGG-ACAACG ' -48.4 ± 



CGCGGAACC/TGTTGC 



-52.3 ± 



3.1 -133.1 ± 10.4 
2.0 -145.9 ± 6.5 



-7.1 ± 0.1 
-7.1 i 0.0 



52.3 
50.9 



-AAGCCTTGT-CCAACG 
^CGCGGAACA/GGTTGC 



-58.0 ± 2.1 -160.2 ± 6.4 -8.4 ± 0.2 57.1 
-50.3 ± 2.1 -135.7 ± 4.7 -8.2 ± 0.6 59.5 



C 



AAGCCTTGC-GCAACG « ^7.0 ± 



CGCGGAACG/CGTTGC 



-48.2 ± 



3.9 -129.0 ± 12.9 -7.0 ±0.1 51.7 
0.9 -132.8 ± 2.9 -7.0 ± 0.0 51.5 



C 



AAGCCTTGG-CCAACG 
CGCGGAACC/GGTTGC 



-57.3 ± 5.4 -156.4 ± 17.0 -8.8 ± 0.3 60.1 
-46.9 ± 2.7 -123.9 ± 5.8 -8.5 ± 0.9 63.6 



C 



AAGCCTTGC-CCAACG 
CGCGGAACG/GGTTGC 



-55.3 ± 
-51.4 ± 



2.6 
1.5 



-152.7 ± 
-140.2 ± 



8.1 
3.4 



-7.9 ± 0.1 
-7.9 ± 0.5 



55.6 
56.7 



^AAGCCTTGG-GCAACG 
^CGCGGAACC/CGTTGC 



c -45.9 ± 
** -52.6 ± 



1.0 -124.4 ± 3.6 .7.3 ± 0.1 54.9 
3.5 -146.0 ± 11.3 -7.4 ± 0.1 52.7 



Imcrfacgs ^th dangling ends' 



C 



AAGCCTTGC-GCAACG 
CGCGGAACG/CGTTGC 
A 



c -45.2 ± 
" -43.7 ± 



2.9 
3.4 



-123.8 
-119.1 



9.4 
11.1 



•6.8 
-6.8 



0.2 
0.1 



50.8 
51.5 



C 



AAGCCTTGC-GCAACG « -46.7 ± 



CGCGGAACG/CGTTGC 
T 



-46.0 ± 



3.7 
3.3 



-128.7 ± 
-126.2 ± 



11.8 
7.5 



-6.8 ± 
-6.8 ± 



0.1 
1.0 



50.9 
51.2 



AAGCCTTGC-GCAACG « -48.6 ± 



CGCGGAACG/CGTTGC 

A 



-47.4 ± 



3.0 
3.2 



-133.7 ± 
-129.8 ± 



9.4 
10.5 



-7.2 ± 
-7.2 ± 



0.1 
0.1 



52.7 
53.1 
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Table S2: Continued. 













> 




3^ 








(kcal / mol) 


(cal / mol K) 


(kcal / mol) 




^AAGCCTTGC-GCAACG 
^CGCGGAACG/CGTTGC 

T 


C 

d 


-46.3 
-51.6 




2.6 
3.7 


-127.9 
-145,4 


± 


8.7 
12.1 


-6.6 
•6.5 




0 

0.1 


43.6 
42.3 


^AAGCCTTGG-TCAACG 
^CGCGGAACC/AGTTGC 

A 


c 
d 


-49.2 
-45.1 




3.5 
2.9 


-136.6 
-123.1 


± 
± 


Il.O 
9.4 


-6.9 
-6.9 


± 


0.2 
0.1 


44.8 
45.7 


^AAGCCTTGC-GCAACG 
^CGCGGAACG/COTTGC 

A A 


c 
d 


-39.7 
-35.5 




4.6 
2.9 


-109.1 
-95.0 


± 


15.2 
5.9 


-5.8 
-6.0 


± 


0.2 
1.1 


38.2 
39.8 


-AAGCCTTGC-GCAACG 
^CGCGGAACG/CGTTGC 

A T 


c 
d 


-44.7 
-43.2 




3,3 
3.7 


-126.0 
-121.0 


± 


ILO 
12.2 


-5.6 
-5.7 


± 
± 


0.2 
0.2 


36.6 
37.2 


^AAGCCTTGC-GCAACG 

T T 


c 
d 


-38.5 
-33.9 




3.0 
1.3 


-105.1 
-89.6 


± 


10.0 
4.5 


-6.0 
-6.1 


± 
± 


0.2 
0.1 


39.3 
41.0 


^AAGCCTTGC-GCAACG 
^CGCGGAACG/CGTT G C 

T A 


c 

A 
O 


-43.3 
-43.9 


± 


9.0 
1.0 


-118.5 
-120.5 


± 
± 


30.1 
3.4 


-6.5 
-6.5 


± 
± 


0.3 
0.0 


43.1 
43.2 


Interfaces with extra centra! nucleotide 
























-AAGC CTTGCAGCAACG 
^CGCGGAACG/CGTTGC 


C 

d 


-43.8 ± 
-45.0 ± 


9.4 
3.5 


-120.2 
-124.1 


± 


31.0 
11.5 


-6.5 
-6.5 




0.2 
0.1 


43.0 
42.7 


^AAGCCTTGTACCAACG 
^CGCGGAACA/GGTTGC 


c 
d 


-45.0 
-45.0 


± 
± 


2.8 
4.0 


-124.5 
-124.8 


± 
± 


9.1 
13.4 


-6.3 
-6.3 


± 


0.2 
0.2 


41.6 
41.4 
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Table S2: Continued. 



Olisomgrs 



T C A A C G 
AGT T G C 



(kcal / mol) 



c -40.6 ± 
" -36.5 ± 



.2.5 
1.9 



AS" AG°3, Ts/ 

(cal / mol K) (kcal / mol) ^--^c i 



■115.4 ± 9.1 
•101.2 ± 6.3 



-4.8 = 0.3 30.5 
-5.1 ± 0.1 31.9 



A C A A C G- 
T G T T Cj C 

GC AAC G 
CGTTGC 

C C A A C G 
GGT T G C 

TGT T G C 
AC A A C G 



-38.3 ± 2.0 -106.8 ± 6.5 

-33.9 ± 1.6 -91.8 ± 5.3 

-43.9 ± 2.4 -121.6 ± 7.4 

-41.1 ± 2.1 -112.6 ± 6.9 

-41.0 ± 2.6 -114.1 ± 8.9 

-36.2 ± 0.9 -98.2 ± 1.8 

-38.8 ± 3.8 -106.6 ± 12.1 

-35.5 ± 2.8 -95.7 ± 9.2 



-5.2 ± 0.2 33.4 

-5.4 ± 0.1 34.6 ' 

-6.2 ±0.1 40.7 

-6.2 ± 0.1 41.2 

-5.6 ± 0.2 36.4 

-5.8 ± 0.3 37.8 

-5.7 ± 0.1 37.5 

-5.8 ± 0.1 38.4 



AGT T G C 

TC A A C G 



-37.1 ± 2.9 -101.2 ± 9.4 
-36.9 ± 2.4 -100.8 ± 8.0 



-5.7 ± 0.1 36.9 
-5.7 ± 0.1 36.9 



* Tm calculated for 4x10"* total strand concentiation. 

The top strand of each system is conventionally represented in the 5' to 3' orientation. Nucleotides 
involved in coaxial stacking interfiaces are represented in bold. 

' Parameters obtained from averaging fits of melting curves. Reported errors are standard deviations 

in the precision of the fitted data. 
Parameters obtained from Tm' vs. \n(Ci/4) plots. Reported errors are standard deviations in the 

precision propagated from the slope and intercept of the 1/Tm vs. in (Cx/4) plot, 
(i). (U), (Hi) ^ ln(CT/4) plots for these systems are shown in Figure S 1 . 
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WHAT IS CLAIMED IS: 



1 1. A method for predicting nucleic acid hybridization 

2 thermodynandcs, the method comprising: 

3 providing a database of thermodynamics parameters; 

4 receiving hybridization information which represents at least one 

5 sequence; 

6 receiving correction data; 

7 receiving a first set of data which represents hybridization conditions; 

8 and 

9 calculating hybridization thermodynamics including net hybridization 

10 thermodynamics based on the hybridization information, the thermodynamic 

1 1 parameters, the correction data and the first set of data. 

1 2. The method as claimed in claim 1 wherein the hybridization 

2 thermodynamics of individual single stranded, bimolecular and higher order 

3 complexes are statistically weighted in a numerical process and the equilibrium 

4 concentration of each species is output. 

1 3 . The method as claimed in claim 2 wherein the correction data 

2 includes folding correction data. 

1 4. The method as clauned in claim 2 wherein the correction data 

2 includes linear correction data. 

1 5 . The method as claimed in claim 1 wherein the thermodynamic 

2 parameters include DNA thermodynamic parameters. 

1 6. The method as claimed in claim 5 wherein the DNA 

2 thermodynamic parameters include dangling end parameters. 

1 7. The method as claimed in claim 5 wherein the DNA 

2 thermodynamic parameters include coaxial stacking parameters. 
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1 8. The method as claimed in claim 5 wherein the DNA 

2 thermodynamic parameters include terminal mismatch parameters. 

1 9. The method as claimed in claim 1 wherein the thermodynamic 

2 parameters include RNA thermodynamic parameters. 

1 10. The method as claimed in claim 1 wherein the thermodynamic 

2 parameters include hybrid DNA/RNA thermodynamic parameters. 

1 11. The method as claimed in claim 1 wherein the thermodynamic 

2 parameters include DNA loop thermodynamic parameters. 

1 12. The method as claimed in claim 1 wherein the hybridization 

2 information represents top and bottom strand sequences which form a duplex and 

3 wherein the hybridization thermodynamics are calculated for the duplex. 

1 13. The method as claimed in claim 1 wherein the hybridization 

2 information represents at least a section of a target and a length of at least one 

3 primer or probe complimentary to the target. 

1 14. The method as claimed in claim 13 wherein the hybridization 

2 thermodynamics are calculated for a plurality of primers or probes complimentary 

3 to the target. 

1 15. The method as claimed in claim 1 wherein the hybridization 

2 mformation represents at least a section of a target and a primer or probe. 

1 16. The method as claimed in claim 15 wherein a length of the 

2 target is longer than a length of the primer or probe and wherein the hybridization 

3 thermodynamics are calculated for a best target/primer or target/probe complex and 

4 for competitive mismatch complexes. 
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1 17. The method as claimed in claim 14 wherein hybridization 

2 information represents at least a section of a target and a primer or probe and 

3 wherein a length of a target is longer than the length of the primer or probe and 

4 wherein the hybridization thermodynamics are calculated for a best target/primer 

5 or target/probe complex and for competitive target/primer or target/probe 

6 complexes. 

1 18. The method as claimed in claim 2 further comprising, 

2 calculating concentration of each species in a solution at a plurality of temperatures. 

1 19. The method as claimed in claim 18 wherein hybridization 

2 information also represents a primer or probe and wherein the length of the target 

3 is longer than a length of the primer or probe and wherein the hybridization 

4 thermodynamics are calculated for a best target/primer or target/probe complex and 

5 for competitive mismatch complexes and wherein the method further comprises 

6 calculating concentration of every species in a solution at a plurality of 

7 temperatures. 

1 20. The method as claimed in claim 19 wherein the hybridization 

2 thermodynamics are calculated for at least two best target/primer or target/probe 

3 complexes and for their corresponding competitive mismatch complexes and 

4 wherein the method further comprises correcting for any interactions between the 

5 at least two best target/primer or target/probe complexes and then: components. 

1 21. A system for predicting nucleic acid hybridization 

2 thermodynamics, the system comprising: 

3 a database of thermodynamics parameters; 

4 means for receiving hybridization information which represents at 

5 least one sequence; 

6 means for receiving correction data; 

7 receiving a first set of data which represents hybridization conditions; 

8 and 
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9 means for calculating hybridization thermodynamics including net 

10 hybridization thermodynamics based on the hybridization information, the 

11 thermodynamic parameters, the correction data and the first set of data. 

1 22. The system as claimed in claim 21 wherein the hybridization 

2 thermodynamics of individual single stranded, bimolecular and higher order 

3 complexes are statistically weighted in a numerical process and the equilibrium 

4 concentration of each species is output. 

1 23 . The system as claimed in claim 22 wherein the correction data 

2 includes folding correction data. 

1 24. The system as claimed in claim 22 wherein the correction data 

2 includes linear correction data. 

1 25. The system as claimed in claim 21 wherein the 

2 thermodynamic parameters include DNA thermodynamic parameters. 

1 26. The system as claimed in claim 25 wherein the DNA 

2 thermodynamic parameters include dangling end parameters. 

1 27. The system as claimed in claim 25 wherein the DNA 

2 thermodynauMC parameters include coaxial stacking parameters. 

1 28. The system as claimed in claim 25 wherein the DNA 

2 thermodynamic parameters include terminal mismatch parameters. 

1 29. The system as claimed in claim 21 wherein the 

2 thermodynamic parameters include RNA thermodynamic parameters. 

1 30. The system as claimed in claim 21 wherein the 

2 thermodynamic parameters include hybrid DNA/RNA thermodynamic parameters. 
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1 31. The system as claimed in claim 21 wherein the 

2 thermodynamic parameters include DNA loop thermodynamic parameters. 

1 32. The system as claimed in claim 21 wherein the hybridization 

2 information represents top and bottom strand sequences which form a duplex and 

3 wherein the hybridization thermodynamics are calculated for the duplex. 

1 33 . The system as claimed in claim 21 wherein the hybridization 

2 information represents at least a section of a target and a length of at least one 

3 primer or probe complimentary to the target. 

1 34. The system as claimed in claim 33 wherein the hybridization 

2 thermodynamics are calculated for a plurality of primers or probes complimentary 

3 to the target. 

1 35. The system as claimed in claim 21 wherein the hybridization 

2 information represents at least a section of a target and a primer or probe. 

1 36. The system as claimed in claim 35 wherein a length of the 

2 target is longer than a length of the primer or probe and wherein the hybridization 

3 thermodynanndcs are calculated for a best target/primer or target/probe complex and 

4 for competitive mismatch complexes. 

1 37. The system as claimed in claim 34 wherein hybridization 

2 information represents at least a section of a target and a primer or probe and 

3 wherein a length of a target is longer than the length of the primer or probe and 

4 wherein the hybridization thermodynamics are calculated for a best target/primer 

5 or target/probe complex and for competitive target/primer or target/probe 

6 complexes. 

1 38. The system as claimed in claim 22 further comprising means 

2 for calculating concentration of each species in a solution at a plurality of 

3 temperatures. 
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1 39. The system as claimed in claim 38 wherein hybridization 

2 information also represents a primer or probe and wherein the length of the target 

3 is longer than a length of the primer or probe and wherein the hybridization 

4 thermodynamics are calculated for a best target/primer or target/probe complex and 

5 for competitive mismatch complexes and wherein the system further comprises 

6 means for calculating concentration of every species in a solution at a plurality of 

7 temperatures. 

1 40. The system as claimed in claim 39 wherein the hybridization 

2 thermodynamics are calculated for at least two best target/primer or target/probe 

3 complexes and for their corresponding competitive mismatch complexes and 

4 wherein the system further comprises means for correcting for any interactions 

5 between the at least two best target/primer or target/probe complexes and their 

6 components. 

1 41. A computer-readable storage medium having stored therein 

2 a database of thermodynamics parameters and a computer program which executes 

3 the steps of: 

4 receiving hybridization information which represents at least one 

5 sequence; 

6 receiving correction data; 

7 receiving a first set of data which represents hybridization conditions; 

8 and 

9 calculating hybridization thermodynamics including net hybridization 

10 thermodynamics based on the hybridization information, the thermodynamic 

11 parameters, the correction data and the first set of data. 

1 42. The storage medium as claimed in claim 41 wherein the 

2 hybridization thermodynamics of individual single stranded, bimolecular and higher 

3 order complexes are statistically weighted in a numerical process and the 

4 equilibrium concentration of each species is output. 
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1 43. The storage medium as claimed in claim 42 wherein the 

2 correction data includes folding correction data. 

1 44. The storage medium as claimed in claim 42 wherein the 

2 correction data includes linear correction data. 

1 45. The storage medium as claimed in claim 41 wherein the 

2 thermodynamic parameters include DNA thermodynamic parameters. 

1 46. The storage medium as claimed in claim 45 wherein the DNA 

2 thermodynamic parameters include dangling end parameters. 

1 47. The storage medium as claimed in claim 45 wherein the DNA 

2 thermodynamic parameters include coaxial stacking parameters. 

1 48 . The storage medium as claimed in claim 41 wherein the DNA 

2 thermodynamic parameters include terminal mismatch parameters. 

1 49. The storage medium as claimed in claim 41 wherein the 

2 thermodynamic parameters include RNA thermodynamic parameters. 

1 50. The storage medium as claimed in claim 41 wherein the 

2 thermodynamic parameters include hybrid DNA/RNA thermodynanwc parameters. 

1 51. The storage medium as claimed in claim 41 wherein the 

2 thermodynamic parameters include DNA loop thermodynamic parameters. 

1 52. The storage medium as claimed in claim 41 wherein the 

2 hybridization information represents top and bottom strand sequences which form 

3 a duplex and wherein the hybridization thermodynamics are calculated for the 

4 duplex. 
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1 53. The storage medium as claimed in claim 41 wherein the 

2 hybridization information represents at least a section of a target and a length of at 

3 least one primer or probe complimentary to the target. 

1 54. The storage medium as claimed in claim 53 wherein the 

2 hybridization thermodynanGiics are calculated for a plurality of primers or probes 

3 complimentary to the target. 

1 55. The storage medium as claimed in claim 41 wherein the 

2 hybridization information represents at least a section of a target and a primer or 

3 probe. 

1 56. The storage medium as claimed in claim 55 wherein a length 

2 of the target is longer than a length of the primer or probe and wherein the 

3 hybridization thermodynamics are calculated for a best target/primer or target/probe 

4 complex and for competitive mismatch complexes. 

1 57. The storage medium as claimed in claim 54 wherein 

2 hybridization information represents at least a section of a target and a primer or 

3 probe and wherein a length of a target is longer than the length of the primer or 

4 probe and wherein the hybridization thermodynamics are calculated for a best 

5 target/primer or target/probe complex and for competitive target/primer or 

6 target/probe complexes. 

1 58. The storage medium as claimed in claim 42 wherein the 

2 program further executes the step of calculating concentration of each species in a 

3 solution at a plurality of temperatures. 

1 59. The storage medium as claimed in claim 58 wherein 

2 hybridization information also represents a primer or probe and wherein the length 

3 of the target is longer than a length of the primer or probe and wherein the 

4 hybridization thermodynamics are calculated for a best target/primer or target/probe 

5 complex and for competitive mismatch complexes and wherein the program 
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6 executes the step of calculating concentration of every species in a solution at a 

7 plurality of temperatures. 

1 60. The storage medium as claimed in claim 59 wherein the 

2 hybridization thermodynamics are calculated for at least two best target/primer or 

3 target/probe complexes and for their corresponding competitive mismatch 

4 complexes and wherein the program executes the step of correcting for any 

5 interactions between the at least two best target/primer or target/probe complexes 

6 and their components. 
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Figure 1 
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Hybridization Information 

Top Strand/Bottom Strand Types 
DNA/DNA p] 

Top Strand Sequence 5-3' 

cccaaaaaaaaaaaccg 



Module 1 



Bottom Strand Sequence 5'-3| Jt| 



Use Complement 
*ggtttttttttttgg* 



JB 



Hybridization Conditions 

[Monovalent cation] ; o.'ios J 



[Mg2+] 

Hybridization 
Temperature 

[Top Strand ] 
[Bottom Strand] 



5e-8 



3e-7 



mol/L 
mol/L 

J mol/L 
J mol/L 



Corrections 

Linear Coirection for Micro Qiips 
(AG° 37(microchip) = a x AG^j^Csolution) + b) 

a= LLI^^^ ^= ^ J 

Top Strand Folding Correction 



kcal/mol AH°= zll'lj kcal/m- 



Bottom Sti-and Folding Correction 
^G^37= J— J kcal/mol AH^>= ^ | kcal/m 



f Piidict Thermodynamics ^ Clear Input 



Figure 2a 
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Duplex sequence 

5^-CCCAAAAAAAAAAACCG-3^ 
3''-*GGTTTTTTTTTTTGG*-5^ 

Experimental conditions 



Module 1 . 



Thennodynamic predictions 



Hybridization type = DNA/DNA 
[Top strand] = 0.5E-07 mol/L 
[Bottom strand] = 0.3E-06 mol/L 

hybridization temperature = 37.0 °C 



Corrections 



In 1.000 MNaCl: 

AH° = -119.3kcal/mol 
AS° = -335.8 eu 
AG° 37 0 = -15.14 kcal/mol 
Tj^ = 52.9°C 



Top strand folding: 

AH° = -37.8 kcal/mol 
AG° 37 Q = -2.10 kcaymol 



In 0.1050 M NaCl and 0.0000 M 
MgC12: 

A H° = -119.3 kcal/mol 

A S° = -348.3 eu 

AG° 37 Q = -1 1.29 kcal/mol 

Tj^ = 42.2 °C 

The net hybridization thermodynamics 
is: 

A G ° 37 Q = -8.74 kcal/mol 
Tj^ = 34.9°C 



Note: 

The net hybridization temperature is the temperature at which the concentration of duplex equals half the 
maximum possible concentration of duplex. 

The net free energy is calculated from the net equilibrium constant at the given temperature- 
Knet=[Duplex]/((Ct-Puplex])*(Cb-[Duplex])), where Puplex] is the concentration of duplex Ct is the 
mitial concentration of top strand, Cb is the initial concentration of bottom strand. 



Figure 2b 
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Hybridization Information 

Target/Primer Types 

Target Sequence 5'-3' 



Module 2 



Primer length 

Number of best primers to be displayed 



J3 



Hybridization Conditions 



[Monovalent cation] 

Hybridization 
Temperature 

[Target] 

[Primer] 



_J mol/L 
I mol/L 



rie^6 1 mol/L 



Corrections 

IJnear Con-ection for Micro Chips 
(AG'^ 3^(microchip) = a X AG*^37 (solution) + b) 



a = 



b = 



f " Pred^^^ : P Clear Inp^^^^ 



Figure 3a 
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Module 2 



Target sequence 



5^-ACCGTTTGTA GTCCGTACGA CACATAACGG TGCATTC 



[Top strand] = O.lE-05 mol/L 
[Bottom strand] = O.lE-05 mol/L 

Hybridization temperature = 37.0 °C 

[Na+] = 1.0000 mol/L 

[Mg2+] = 0.0000 mol/L 

The 2 best primers of length 15 are: 



From position 28 to 42: 5'-GGTTGCAATGCACCG -Z' 

A H° = -132.0 kcal/mol AS° = -355.6 eu AG° 37 g = -21.71 kcaymol Tj^ = 
70.2 °C 

From position 35 to 49: 5^-GCAGCATGGTTGC.Ai\ -3^ 

AH° = -124.8kcal/mol AS° = -336.5 eu AG° 37 q = -20.42 kcaVmol Tj^ = 
68.4 °C 



Experimental conditions 



Corrections 



Hybridization type = DNA/DNA 



No corrections 



Figure 3b 
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Module 3 



Hybridization Information 

Target/Primer type 

Target 5-3' Primer Sequence 5-3' pi 

• acgcttgaatgcagttaatgcc ^ tgaatgckgt " ' 



Minimum percent stability of alternative binding sites , ^- 



compared to the most stable binding site 
Number of base pairs required to compute the solution i^^Z^^^lf 



Hybridization Conditions 




User defined values for [Na+J and [Mga+J p] 


[Monovalent cation] 




mol/L 






mol/L 


Hybridization 
Temperature 

[Target] 

[Primer] 


le- 6 * ■ 


1 mol/L 
1 mol/L 



Corrections 

Linear Correction for Micro CJiips 
(AG^ 37(microchip) = a x AG^g^Csolution) + b) 



a = 



b = 



Submit I Clear In^^ 



Figure 4a 
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Module 3 



Target sequence 



5^-ACGCTTGAAT GCAGTTAATG CC-3^ 



Primer sequence 



3'-TGACGTAAGT-5^ 



Experimental conditions 



Corrections 



Hybridization type = DNA/DNA 



No corrections 



[Top strand] = O.lE-05 mol/L 
[Bottom strand] = O.lE-05 mol/L 

Hybridization temperature = 37.0 °C 

[Na+] = 1.0000 mol/L 

[Mg2+] = 0.0000 mol/L 

Number of base pairs required to compute the 

solution = 5 

Best primer site 



from target position 8 to position 17 

GAATGCAGTTAA 
TGACGTAAGT 



A H° = -26.2 kcal/mol AS° = -70.7 eu AG° .70 = -4-28 kcal/mol 



= -9.9 °C 



Figure 4b 
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Hybridization Informatinn 

Target/Primer Types 
^ DN/^DNA JV| 

Target Sequence 5'-3' 



Module 5 



B 



Find best primer in sequence section ranging from 
nucleotide number: 



Primer length 
PI Number of best primer 



to 

15 f 



Percent stabilitty of altemative binding sites compaied to the most stable binding site 50 I 
Number of base pairs required to compute the solution 



Hybridization Conditions 

■ User defined values fij^^^NQj^^^^ 

[Monovalent cation] f i 1 mol/L 

tMg^"^ fo"" I mol/L 
Hybridization ;3-... 
Tempemture ^= ^ 

[Target] rj^ls " • •| moI/L 

[Primer] ^ ie~6 ' " " | moI/L 

I Suijmit 



Corrections 

Linear Coixection foj- Micro Chips 
(AG° 37(microchip) = a x AG^37(solution) + b) 

^i—J b= ,0,:,: j 



a = 



I Clear Input 



Figure 5a 
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Target sequence 

Module 5 



5/- 

AGGTCCATGCTTTGGAACAGCTACTTGAACCGATCATGGACACTGACGGATAAC( 

Experimental conditions Corrections 



Hybridization type = DNA/DNA 

[Top strand] = O.lE-05 mol/L No corrections 

[Bottom strand] = O.lE-05 mol/L 

Hybridization temperature = 37.0 °C 

[Na+] = 1.0000 mol/L 

[Mg2+] = 0.0000 mol/L 

Number of base pairs required to compute the 

solution = 7 

Best primer search area from position 1 to 
position 60 



Best primer # 1: 



from target position 35 to 49 

5 ' -TCATGGACACTGACGGA- 3 ^ 
3 / -GTACCTGTGACTGCC - 5 

AH°=-123.1kcal/mol AS° = -331.5 eu AG° 37 q = -20.27 kcaymol Tj^ = 
68.4 °C 



Best primer # 2: 



from target position IS to 32 
5 ^-ACAGCTACTTGAACCGA- 3 ^ 
3 -GTCGATGAACTTGGC - 5 ^ 



'A H° = - 1 25.0 kcaymol AS° 
66.1 °C 



-339.6 eu AG° 37 q = -19.67 kcal/mol Tj^ = 



Figure 5b 
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initial guess: 



Total 



temperature T==0 C 



Total 




Multiplex PCR Design 



^ Gene I Gene II ^"^^ 
— ^ Gene l ' ' I I M M 



P2 P4 



Figure 7 
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Prediction of Molecular Beacon 
Hybridization 



Duplex Formation 

CCCAAAAAXAAAAACCG 
+ 

GGTTTTTYTTTTTGG 



X-Y 



CCCAAAAAXAAAAACCG 
GGTTTTTYTTTTTGG 



A 

/A 



A-T -13.24 47.9 

A-A -9.94 39.1 

A-C -9.03 36.0 

A-G -10.36 40.4 

A X A 



Beacon Folding 



Random Coil Beacon 



A 
A 
A 

C 



A 
A 
A 
A 
C 



AG»37 =-2.1 



kcal/mol 



Tm = 55.2 °C 



Net Hybridization 

Hairpin Beacon 
+ 

Target 



C C 
C-G 

T-A 

C-G 

G-C 

C-G 

4i 

Hairpin Beacon 



1/ 



-CGTCCCCAAAAAAAAAAACCGACG-C 
\ 3 ' -GGTTTTTXTTTTTGG- 5 ' 



Target 



X-Y 



AGO37 (Effective) 



(Effective) 



3'GG I I I I I I I I I I I GG^' 
3'GGI I I I lAI I I I IGG^' 
3'GGI I I I ICI I I I IGG5' 
3'GGI I I I IGI I I I IGG^' 



A-T 
A-A 
A-C 
A-G 



Exp. 
-10.49 
-6.66 
-6.72 
-7.62 



Pred. 
-10.69 
-7.39 
-6.48 
-7.81 



Exp. 
42 
27 
23 
28 



Pred. 

42.4 

26.8 

21.1 

29.5 



0.105 MNaCl 0.001 MMgClj [beacon] = SxlO-^M [target] = 3x1 0-''M 
Bonnet et al. (1999), Proc, Nat. Acad. Set USA 96, 6171-6176 



Figure 8 
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Simulation of Molecular Beacon Hybridization 




♦ F(RC) 

■ F(Dup) 

A F(Hairpin) 

X F(target bound) 



120 



Temperature 



Figure 9 
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Match VS. Mismatch Hybridization 



Match Site 



Mismatch Site 



Probe 



Probe^M 



Target DNA 



Equilibria 



Match Hybrid T + P 



T-P 



Mismatch Hybrid T 4- P^ 



Double Hybrid T + 2P 



'p.pMM 



X-P-P^^M 



[T-P] 



M [-pj 

[■p.p.pMMj 

Kdh= [T] [P]2 



' Given CTargetCtotal], Cprobe[total]and the 3 equilibrium constants above, 
it is trivial to solve for the concentrations of all species 

• Since AC^j and AH° are known, calculate K's at all temperatures 

• Simulate hybridization at all temperatures - optimize specificity 

• More complex model would also include single-strand folding equilibria 

FigurelO 
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Match vs. Mismatch Hybridization Simulation 



c 
o 

(0 



o 
o 
c 
o 
o 

i 

(0 

E 
1- 

o 

z 



1.00 
0.90 
0.80 
0.70 
0.60 
0.50 
0.40 
0.30 
0.20 
0.10 
0.00 




♦ Fraction Target in R.C. 
■ Fraction Probe in R.C. 
A Fraction Match Hybrid 
O Fraction Mismatch Hybrid 
▲ Fraction Double bound 




20 



40 



60 



80 



100 



Temperature (^C) 



Figurell 
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